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ABSTRACT 



In quantum information theory, a geometric approach, known as "quantum information 
geometry," has been considered as a powerful method. In this thesis, we give a computa- 
tional geometric interpretation to the geometric structure of a quantum system. Especially 
we introduce the concept of the Voronoi diagram and the smallest enclosing ball problem 
to the space of quantum states. With those tools in computational geometry, we analyze 
the adjacency structure of a point set in the quantum state space. Additionally, as an 
application, wc show an effective method to compute the capacity of a quantum channel. 

In the first part of this thesis, wc show some coincidences of Voronoi diagrams in a 
quantum state space with respect to some distances. That helps us to reinterpret the 
structure of the space of quantum pure states as a subspace of the whole space. More 
properly, we investigate the Voronoi diagrams with respect to the divergence, Fubini-Study 
distance, Bures distance, geodesic distance and Euclidean distance. 

For one qubit (or two level) states, the whole space is expressed as the Bloch ball. 
In the Bloch ball, we analyze the coincidence of the Voronoi diagrams in two different 
settings: l)the diagram in pure states when Voronoi sites are taken as pure states and 
2)the diagram in mixed states when sites are taken as pure states. We show that in both 
cases, all the diagrams coincide. This clear result is because of the symmetry specific for 
one qubit states. 

For three or higher level systems, we investigate the diagrams in pure states. The nat- 
ural embedding of the quantum state space into a Euclidean space is no longer symmetric 
as in one qubit case. Consequently the coincidence of the Euclidean Voronoi diagram 
and the divergence- Voronoi diagram does not hold in a higher level system. However the 
coincidence of the divergence-, Fubini-Study- and Bures- Voronoi diagrams still holds. 

In the second part, we propose a method to compute a capacity of a quantum communi- 
cation channel and show the result of the actual computation. Wc show that our method is 
sufficiently effective not only for one qubit states but for three level states. It is a practical 
application of the theoretical result shown in the first part; the theorems in the first part 
guarantee the correctness of the algorithm used in the second part. The algorithm uses 
Welzl's algorithm to compute a smallest enclosing ball. Although the original algorithm 
introduced by Welzl is only for the Euclidean space, we show that the same method is 
useful for non-Euclidean space. We also implement the algorithm and experiment it to 
prove it is practical. 
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Chapter 1 



Introduction 



In this chapter, we explain the background for our research and the summary of our 
contribution. Our contribution has mainly two aspects: computational geometry 
and quantum information theory. 

In computational geometry, our contribution is shortly described as introduc- 
tion of another Voronoi diagram with a distortion measure and an algorithm to 
solve the smallest enclosing ball problem in that measure. 

In quantum information theory, we contribute to re-interpret the structure of 
a quantum state space to some extent, and proposed a practical algorithm to 
compute the capacity of a quantum channel. 

We explain summarized background for topics related to our contribution. The 
explanation of the background is divided into two parts: one is about computational 
geometry and the other is about quantum information theory. Then, we explain 
the outline of the contribution of this dissertation, and show how this dissertation 
is organized. 

1.1 Computational geometry 
1.1.1 Voronoi diagrams 

Voronoi diagrams and Delaunay triangulations have had an important role in com- 
putational geometry. Voronoi diagrams are not only useful for such applications 
as numerical calculation and visualization, but also useful for theoretical interpre- 
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tation of a geometric object. 

Practically, Voronoi diagrams are used in many fields such as geophysics, me- 
teorology, astrology, geometric information system, city planning, and so on. Ad- 
ditionally, nowadays one of the most important applications of Voronoi diagrams 
is computer vision. In particular, because of the recent rapid development of com- 
puter graphics in entertainment media such as movies and games, Voronoi diagrams 
are getting more and more important. Other practical examples of applications 
are explained in [87j . 

A Voronoi diagram is a division of a space. The applications listed above all 
deal with a geometric space and need to divide it so that it can be computed in 
reasonable time. The essence of Voronoi diagrams can be explained as follows. 
Suppose that some points (which we call sites) are given and you want to divide 
the space into some regions so that each region expresses dominance of a point. 
If for any point in a region, the nearest site is the site included in the region in 
problem, then the division is called a Voronoi diagram. A Delaunay triangulation 
is a dual of a Voronoi diagram; for any two sites, draw a line between them if there 
is a Voronoi edge between them, and you obtain a Delaunay triangulation. Fig. 11.11 
shows a example of Voronoi diagram and Delaunay triangulation. 

As Voronoi diagrams used in the wider area, the more complicated diagrams 
with many Voronoi cells became needed and consequently, the robustness of the 
computation became emphasized. Actually, the general real world problems which 
needs Voronoi diagrams are badly positioned so that the resulting diagrams are 
degenerating. Degeneracy of Voronoi diagrams happens when the sites are cocir- 
cular. The evilness of degeneracy lies on the weakness for perturbation; only slight 
move of a site will change the topological position (See Fig. II. 2p . In calculation 
in a computer where numbers are not expressed rigidly, this kind of evilness may 
cause inconsistent situation such as "A point is geometrically in a certain area, but 
topologically out of the area." One off the researches to overcome this problem is 
by Sugihara and Iri |112^ I113j . They achieved the robustness by focusing only on 
topologies and ignoring the geometries when computing about the relative position 
of edges. Thanks to those researches about the robustness of actual computation, 
now Voronoi diagram with billions of sites is shown to be computed. Isenburg et 
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Figure 1.1: An example of Voronoi diagram and Delaunay triangulation: Solid 
lines are a Voronoi diagram, and dotted lines are a Delaunay triangulation 

al. [58] showed an algorithm to calculate a tessellation of terrain data with a billion 
of triangles using some heuristics. 

The significant application which emerged recently is fine art. Fritzsche et 
al. proposed a algorithm to synthesize an authentic look mosaic structure from 
a picture. Sugihara |111| synthesize a artistic pattern. Since there is often a fractal 
structure behind a beauty in the nature, he proposed an algorithm to create a 
nature-like shape by combining Voronoi diagrams with fractal. 

1.1.2 Generalized Voronoi diagram 

In the original "normal" Voronoi diagrams, the sites are given as points and Eu- 
clidean distance is used to decide the dominance of each region. The generalization 
of Voronoi diagrams mainly goes into two ways: a) define a site as a set of point 
instead of one point, or b) use a general distance instead of the Euclidean distance. 
In the direction to consider a general site, Voronoi diagrams for line segments 
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Figure 1.2: An example of degenerate Voronoi diagram: Sites are in a cocircular 
position and Voronoi edges meet at one point. 

were intensively researched from the late 70's, by Drysdale and Lee ^9], Drys- 
dale [28], Kirkpatrick [72], Lee and Drysdale |;73j, Imai et al. |!57j, Sharir |106j . 
Fortune [35], Yap [124j . Clarkson and Shor [20], Goodrich et al. [38], Burnikel et 
al. [T7], Rajasekaran and Ramaswami [981 199j . and Deng and Zhu [25]. The rather 
recent ones are Voronoi diagrams for circles or balls. Concerning the Voronoi di- 
agrams whose sites are given as circles, Kim et al. |641 165] proposed a algorithm 
which is robust even for the degenerating case. 

Here, note that when we say "distance" in the context of Voronoi diagrams, it 
does not necessarily satisfy the axioms of a distance. In this dissertation, we use 
the word "measure", or "pseudo-distance" in a confusing context. 

Another direction is to use a general distance in a general space. One of the 
simplest in this direction is Voronoi diagrams for the weighted distance [87|. It 
means that each site has a weight and the distance is measured according to the 
weight. If Wi is the weight for the site Sj, then the multiplicatively weighted dis- 
tance to the point x is defined as dweighted(si) s;) = \x — Si\ /wi. Fig. 11.31 shows 
an example of multiplicatively weighted Voronoi diagrams. The Voronoi edges for 
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Figure 1.3: An example of weighted Voronoi diagram: Each number assosiated to 
the sites means its weight 

a multiphcatively weighted Voronoi diagram becomes a part of a circle, so called 
ApoUonius circle. 

The weighted distance is distorted to some extent, but it still based on Euclidean 
distance. Using more general distance, more distorted or sometimes pathological 
Voronoi diagrams can be obtained. For example, Onishi and Itoh investigated 
Riemannian Voronoi diagram [91]. In classical information theory, Onishi and 
Imai [iU [90] and Nielsen et al. [80] are for divergences. The detail of divergence 
Voronoi diagrams are explained in the next section. 

The combination of the two directions described above can also be considered. 
Generalizing the way to decide the polygon mesh used in computer vision, Asano 
introduced an aspect-ratio Voronoi diagram and analyzed its computational com- 
plexity and Asano et al. also introduced an angular Voronoi diagram [7J. In 
those diagrams, Voronoi sites are given as line segments and the distance is given 
as a visual angle or the aspect ratio of the triangle composed by a point and a line 
segment respectively. In both cases, the Voronoi edges are curves of degree three 
and the regions can be complicated; the region dominated by the same line segment 
can separated by points. Fig. 11.41 shows the example of angular Voronoi diagrams. 
Asano et al. also generalized those Voronoi diagrams to obtain an abstract notation 
of Voronoi diagrams and analyzed the characteristics of those diagrams [8] . 

When we think about actual computation of those generalized distortion 
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Figure 1.4: An example of an angular Voronoi diagram (drawn by H. Muta): 
Although it can be drawn in two dimensional space, its measure is distorted and 
the Voronoi edges are generally cubic curves 

Voronoi diagrams, the arising problem is how we can achieve its robustness. Muta 
and the author analyzed the extended the notion of degeneracy of Voronoi diagrams 
from the viewpoint of the computational robustness [79]. The main discussion in 
the paper is about the conditions to make the number of crossing points of Voronoi 
edges jump with perturbation. It is only about angular Voronoi diagrams, but the 
concept of degeneracy in the paper is also applicable to general Voronoi diagrams, 
although there is only few attempt to actually compute the diagrams explicitly. 

There are still other Voronoi diagrams with some "strange" measure spaces. In 
the book by Okabe et al. [,87j, many of them are introduced with analysis about 
their computational complexity. 

1.1.3 Voronoi diagrams in classical and quantum information 

Information theory is considered to have been founded by Shannon |105j . He 
showed the bound for the capacity of a channel by coding the source message 
of the channel. Its coding strategy is decided by the probabilistic distribution 
of the source message. Thus, information theory is mainly based on probability 
theory and statistics. The most important quantities used in this field are entropy. 
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and Kullback Leibler divergence (or relative entroby). Kullback-Leibler divergence 
is defined as some kind of "distance" of two probabilistic distribution. Hence, 
geometry in information space can be considered but its structure is very distorted 
and far from intuition. Kullback-Leibler divergence does not satisfy the axioms of 
distance (so it is not a distance in a rigid sense); for example it does not satisfy 
the law of triangle inequity, or is not commutative either. However, its distorted 
and strange properties give a rich field of a research for computational geometry. 

A computational geometric analysis was done by Onishi and Imai |89[ |2Dj , On- 
ishi [88], and Sadakane et al. |101j . A Voronoi diagram and Delaunay triangulation 
are defined with respect to the Kullback-Leibler divergence, and are shown to be 
the extensions of the Euclidean counterparts. The Voronoi diagram is computed 
from an associated potential function instead of a paraboloid which is used in a 
Euclidean Voronoi diagram. 

In the same line, Nielsen et al. [80] showed some properties of Voronoi diagram 
with respect to Bregman divergence, which is generalization of Kullback-Leibler 
divergence. Using the Voronoi diagram, Nielsen et al. [ST] also showed that Welzl's 
algorithm to solve the smallest enclosing ball problem is also applicable to Bregman 
divergence. 

We extend the Voronoi diagram in classical information to the quantum world. 
In quantum information theory, there is a natural extension of the Kullback-Leibler 
divergence, and it is called a quantum divergence. We introduce a Voronoi diagram 
with respect to the quantum divergence, and analyze its structure. Additionally 
we consider other diagrams with respect to some distances. Comparing the dia- 
grams, we can compare the structures of some distance spaces and consequently 
some problem concerning a certain distance can be replaced by another problem 
of another distance. 

1.1.4 Smallest enclosing ball problem 

The smallest enclosing ball is namely a problem to compute the smallest ball which 
contains given points. It has variety of applications; collision detection, facility 
location, automated manufacturing, and so on. It is a geometric problem but has 
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some aspect of combinatorial optimization. 

The first theoretically effective algorithm was given by Megiddo [77|. In spite 
of its astonishing idea of pruning- and- search, Megiddo's algorithm was impracti- 
cal because there is a big constant hidden behind a big-0 notation. Welzl jl21j 
gave the first practical algorithm based on Seidel's randomized linear programming 
algorithm jl04j . In these algorithms, however, the complexity is the exponential 
of the dimension. Eventually Matousek et al. [75] discovered subexponential time 
algorithm. The most efficient algorithm known so far is Fischer and Gartner's algo- 
rithm [34j. They gave an 0(d^(1.438)'^)-time algorithm. Fischer also implemented 
a program to compute the smallest enclosing ball problem as a part of CGAL |19j . 

Nishitoba et al. [83] connected this line of the research to the computation 
of the Holevo capacity in quantum information theory. He also shed light on 
combinatorial aspect of this problem. This direction is followed by Nishitoba [82j to 
analyze the combinatorial structure. Although Hayashi et al. had already showed 
the method using the smallest enclosing ball problem, they first mentioned the 
necessity of the fast algorithm with respect to the dimension to extend the existing 
method to the higher level system. Actually the dimension of the space when we 
think of the smallest enclosing ball problem is — 1 for d-level system; this grows 
too rapidly from the viewpoint of practical computation. 

1.2 Quantum information theory 

1.2.1 Quantum computation and quantum information 

Feynman [33] is considered to be one of the earliest to show an idea to apply quan- 
tum mechanics to computation. His idea comes from the fact that in quantum 
mechanics, huge amount of computation is needed to compute the behavior of par- 
ticles; he considered that particles which behave according to the law of quantum 
mechanics can be used to compute a quantum behavior itself. 

Deutsch [26] is the first to show this idea is really useful in some problem. 
He showed a problem which can be solved by a quantum computer exponentially 
faster than a classical computer. Although the problem proposed by him is rather 
artificial and it is only to show the computational gap between a classical computer 



8 



and a quantum computer, it has a significant meaning as a first example to show the 
power of superposition of quantum states in computation. The original algorithm 
proposed by Deutsch is only for one bit but it was shown to be extended to n-bit 
by Deutsch and Jozsa [27J, and Cleve et al. [2l] gave another improvement . 

One of the greatest works in quantum computation was by Shor [llOj. He 
showed that using a quantum computer, factoring of integer and discrete loga- 
rithms can be solved in a polynomial time. It is the first example in which quan- 
tum computer is exponentially faster than classical computer. It was sensational 
because the difficulty of factoring is a guarantee for the security of the existing 
public key cryptosystem. 

Another famous algorithm for quantum computer is database search algorithm 
by Grover [39j. His algorithm is quadratically faster than classical one. His orig- 
inal algorithm is improved and generalized by Grover himself |4H HO] and Biham 
et al. [13]. Another algorithm is about integration. Abrams and Williams [1] 
showed an algorithm for multi-dimensional integration. Analysis on some classes 
of functions is done by Novak \8^ and Heinrich [39] . For quantum algorithms for 
numerical integration, surveys are written by Heinrich [481 150j . 

Miyake and Wadati [78] showed that the Fubini-Study distance in a quantum 
state space has a special meaning for quantum search algorithm. In the continuous 
version of Grover algorithm, a quantum state follows a path which is geodesic in 
the Fubini-Study distance. 

Quantum information theory has been considered as a primitive backbone for 
the quantum computation. The typical theme of quantum information is "What 
can be possible using a quantum channel?" Its difficulty is based on the charac- 
teristics of the measuring of quantum states. Even for the two completely same 
quantum states, the result of the measurement may be different and may distribute 
probabilistically. Consequently, the typical objective of quantum information the- 
ory is to distinguish some different quantum states by measuring. 

For quantum information theory, the invention of the quantum cryptosystem 
[11] is an important epoch-making event; it became a trigger to make active the 
research in that field although it was not the very start of quantum information 
theory. From the practical point of view, quantum cryptosystem is believed to be a 
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very near to utilization in the real world. Tajima et al. |114| reported a success of an 
experiment toward realization of quantum cryptosystem with a realistic settings. 
Additionally, some venture companies, such as MagiQ Technologies ^74j and id 
Quantique [56], are emerging in this field. 

Some aspect of quantum information theory is to investigate a kind of dis- 
tance between two different quantum states. Depending on the situation, several 
distances are defined in quantum states. In quantum information geometry, the 
structure of those distances is researched |3j 95\ . 

The quantum divergence have been used as an informational distance from a 
quantum state to another. In particular, it played an important role in quantum 
hypothesis testing |5H I85j and an estimation of a capacity of a quantum channel [52\ 
[53] . This informational measure is the main of our interest. It is not symmetric and 
so distorted that its structure is difficult to understand. We have been motivated 
to understand its structure and clarify its geometric properties. 

1.2.2 Power of entanglement and additivity conjecture 

Entanglement is considered to be one of the most important and interesting objects 
in quantum information theory, and actually provides a hot field of research. As is 
described above, the result of the measurement of a quantum state may distribute 
probabilistically; let us compare it to the coin tossing. Then entangled states are 
like correlated coins; the probability whether one will show the top or tail is related 
to the result of the toss of the others. 

This correlation, a strange behavior of particles, were pointed out by Einstein 
et al. [31] and the claim in this paper is known as Einstein-Podolsky-Rosen (EPR) 
paradox. Its original intention was to show the paradox of quantum mechanics. If 
there were such correlation, the locationally separated particles can provide a mean 
to convey information more rapidly than the speed of the light; they claimed it is 
a contradiction. In spite of its original intention, their result is now known to be 
a fundamental principle of quantum mechanics. Since the correlation given by the 
EPR paradox was known to be a break of the inequality shown by Bell which shows 
the necessary condition for a given artificially made probabilistic distribution to 
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be really realizable. Some of the recent research concerning Bell's inequity is by 
Tsirelson [IISl UHl EZ] , Avis et al. [TO], and Ito [59] 

One of the lines of the research concerning entanglement is based on the rather 
philosophical question: "What is communication?" It is shown that sharing entan- 
gled states helps two parties to win some sort of games. The researches in quantum 
games are based on the idea that whether ones are communicating or not is judged 
by whether they can do something they could not do without any share of infor- 
mation. Avis et al. f9j showed that two parties sharing entangled states can make 
a good performance in the graph coloring game compared to the case that there 
is no communication. This result indicates that they are surely communicating 
something, although it is much weaker than classical communication. It is called 
a pseudo-telepathy because it looks like a telepathy but can win in only limited 
games. 

Another direction of the research is to evaluate numerically how much the 
states are entangled. One direction is to measure some kind of distance from the 
maximally entangled state to the state in problem. Some entanglement measures 
were proposed by Bennett et al. [12]. The generalization of the measures which 
means investigation for the condition which the entanglement measures must satisfy 
is done by Vidral et al. [HQl [JTOl EH], and Rains [97]. 

"How much do entangled states contribute to the capacity of a quantum com- 
munication channel?" has been considered as an important problem. The problem 
can be described more precisely as follows: does the Holevo capacity of a given 
channel make any difference depending on whether its domain is restricted to sep- 
arable states (i.e. not entangled states) or not? In a mathematical sense, such a 
problem is stated as an "additivity problem." It is conjectured that the additivity 
holds for any quantum state space. In other words, it is believed that entangled 
states give no power to the quantum channel with respect to some measures. 

Concerning the problem of sending a classical message via a quantum chan- 
nel, Holevo showed the upper bound for its capacity [52i 153] . Holevo [5l] and 
Schumacher-Westmoreland [102] independently showed theoretically that the up- 
per bound can be attained. 

Shor [109] proved that some open problems concerning the additivity with re- 
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spect to some measures are all equivalent. In particular, they proved the equiva- 
lence of the additivity of the Holevo capacity and the additivity of the minimum 
entropy output. The additivity of the minimum entropy output is equivalent to 
the limit of the multiplicativity of the p-norm as p — > 1. 

Although the conjectures are not solved completely, it is confirmed to hold for 
some classes of channels. About unital channels. King [67j proved for unital qubit 
channels; Fujiwara and Hashizume [37j, King ^9], and Amosov [4j for depolariz- 
ing channel; Matsumoto |76j . Datta et al. [23], and Alicki [2] for Werner-Holevo 
channels; Fannes et al. [32] and Datta et al. [22] for the transpose depolarizing 
channel; Datta and Ruskai [24j for some asymmetric unital channel. About non- 
unital channels, Shor [108] and King ^68] proved for entanglement-breaking chan- 
nels; Wolf [123] for a modification of the Werner-Holevo channel; and King [70] for 
diagonal channels. 

Although it is also conjectured that the multiplicativity of the p-norm holds, a 
counterexample for p > 4.79 is discovered by Werner and Holevo [122j . However, 
it is still believed it holds for a sufficiently small p. King and Ruskai [71] showed 
a condition under which the multiplicativity holds if p = 2. The condition shown 
there holds for typical examples, and accordingly, it is a strong support for the 
multiplicativity conjecture especially when p = 2. On the other hand, almost 
nothing is known around p = 1, and thus, the additivity problem is considered to 
be extremely difficult. 

Those problems related to entanglement also motivated us. Our computational 
geometric approach clarifies the structure of an entanglement measure. Our nu- 
merical computation for the Holevo capacity is also related. If there were a time 
effective and numerically robust algorithm to compute the Holevo capacity even 
for a high level system, it would be a good tool to check the additivity conjecture. 
Our algorithm is a step toward it. 

1.2.3 Numerical estimation of a quantum channel 

Generally, a space of quantum states has a complicated structure. For d-level 
system, the whole space is known to be closed convex object in a Euclidean space of 
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real dimension d? — l. Kimura [66] showed the hst of the inequahties which formahze 
the conditions for embedding of a quantum state space into EucHdean space and 
it is convincing the comphcatedness of the space. Because of the comphcatedness, 
the capacity of a quantum channel is difficult to compute although it is important 
for quantum related engineering. Whether there is a practically efficient algorithm 
to compute the capacity is one of our interests. 

Hayashi et al. [47] and Oto et al. [93^ [94] showed an effective method to numer- 
ically compute the Holevo capacity of one-qubit quantum channel. With an actual 
numerical computation, Hayashi et al. [47j showed that there is a case that needs 
maximal number of points to determine the smallest enclosing ball; this means 
the quantum state space with respect to the divergence is distorted compared 
to Euclidean space. Actually Fig. 11.51 shows an example of a divergence-sphere 
and we can observe it is really distorted. About the three level system, Osawa 
and Nagaoka [92] were the first to show an example of numerical computation. 
They proposed a quantum version of Arimoto-Blahut algorithm |5l [TJ] to numeri- 
cally compute a Holevo capacity, and confirmed that the additivity holds for some 
three-level examples. Our motivation partially comes from their work. If there is a 
faster algorithm to compute the Holevo capacity, although this direction will never 
help us to prove it, we can be more convinced with the additivity conjecture or 
otherwise, can find a counterexample. 

The Holevo capacity is defined as the capacity of a quantum channel when it 
sends classical message. Consider the setting that you send a classical message 
via a quantum channel and suppose that a probabilistic distribution of source 
messages and the way of encoding a message are varying parameter. The Holevo 
capacity is the maximum of conveyed information with those varying parameters. 
Its concept is illustrated in Fig. 11.61 An upper bound for the capacity is proved by 
Holevo [521 [53], and it is proved to be attained [M ] 1102] . 

The method we introduce in this dissertation is based on the algorithm to 
compute the smallest enclosing ball. This is an extension of the method used by 
Hayashi et al. [17] and Oto et al. |93[ . The smallest enclosing ball problem itself 
is also important in computational geometry and is still to be solved from some 
aspect; even in the Euclidean distance space, a polynomial time algorithm in the 
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(a) (b) 
Figure 1.5: An example of a divergence-sphere : 3D view of a divergence-sphere 
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Figure 1.6: An explanation of the setting for the Holevo capacity 



dimension of the space is not known. 

1.3 Contribution of this dissertation 

Our main contribution is a computational geometric interpretation of a quantum 
state space. In particular, wc introduce a concept of Voronoi diagrams and the 
smallest enclosing ball problem in a quantum state space. Those standard tools in 
computational geometry help to clariiy the adjacency structure of a point set in a 
quantum state space. Another aspect of our interpretation is that we show some 
relation between quantum information and quantum computation. Actually we 
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Figure 1.7: A description of a part of our contribution 



show some different distances used in quantum information and quantum compu- 
tation give a same Voronoi diagram. The setting for the problems described above 
and the variety of the (pseudo-)distances in a quantum state space are illustrated 
in Fig. 11.71 Moreover as an application of such a geometric interpretation, we pro- 
pose an algorithm to compute a capacity of a quantum channel, and show, by an 
experiment, that it is really practical. 

By proving coincidences of Voronoi diagrams, we show a connection between 
some distances which were considered differently. An especially significant point is 
that we showed a coincidence of Voronoi diagrams with respect to the divergence 
and the Fubini-Study distance. The divergence is an important measure in quan- 
tum information theory, and is used for quantum hypothesis testing [5H [85] and 
estimation of the capacity of a quantum channel \52\ [53], while the Fubini-Study 
distance gives a convergence path in Grover's search algorithm [78|. We bridge 
those topics which had not seemed to be related, but had been considered to be 
both important for quantum-related researches. Moreover, we also show a connec- 
tion between the Bures distance and the divergence. The Bures distance is used as 
a measure of entanglement |102[ [46]. Although for pure states, the Bures distance 
is fundamentally the same thing as the Fubini-Study distance, it is meaningful to 
know a connection between the Bures distance and the divergence which are used 
in different contexts. 

Another interesting point is that the Fubini-Study distance is only for pure 



15 



states while the divergence is not defined for pure states. To find the relation 
between those exclusive measures, we introduce a natural definition of Voronoi 
diagram with respect to the divergence in a space of pure states. Here, note 
that although the divergence for pure states is not defined, the divergence- Voronoi 
diagram is naturally extended to a space of pure states by taking a topological 
closure. A space of pure states has a simple structure and the Fubini-Study distance 
is defined as a very natural distance in it. It is meaningful because the whole space 
including pure and mixed states has a complicated structure and the quantum 
divergence is a distortion measure in it. Thus we give a connection between a 
simple natural structure and a distorted structure. 

From the viewpoint of computational geometry, our contribution is that we 
introduce Voronoi diagrams in a distortion space and characterize it. Our main 
interest is the quantum divergence — the most distorted one among the measures 
defined in quantum state space. A Voronoi diagram with respect to the quantum 
divergence is a natural extension of a diagram with respect to Kullback-Leibler 
divergence in classical information theory, and we reveal its geometric properties. 

For pure states in the space of one-qubit quantum states, we show the coin- 
cidence of Voronoi diagrams with respect to some distances — the divergence, 
the Fubini-Study distance, the Bures distance, the geodesic distance and the Eu- 
clidean distance [GOlllS]. As an application of this fact, we introduce a method to 
compute numerically the Holevo capacity of a quantum channel [93^ [Ml [47] . The 
effectiveness of this method is partially based on the coincidence of the diagrams. 
Moreover, also as to the diagrams in mixed states, we found the coincidence of 
some of them. The diagrams with respect to the three distances — the divergence, 
the Fubini-Study distance, and the Bures distance — coincide [62l [63] . 

A natural question that arises after this story is "What happens in a higher level 
system?" For a higher level system, the diagrams with respect to the divergence and 
the Euclidean distance do not coincide [61[ 163]. On the other hand, the diagrams 
with respect to the divergence, the Bures distance and the Fubini-Study distance 
still coincide for a higher level. 

We also show that Welzl's algorithm is also applicable to quantum state space. 
Most of the idea of its proof is by Nielsen et al. [80], but they did not mention 
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about quantum divergence. Following their idea, we show that the smallest en- 
closing problem in quantum state space with respect to quantum divergence obey 
to the axioms of LP-type problem which is essential condition to show the Welzl's 
algorithm is effective. 

As an application of the theoretical result we proved, we propose an algorithm 
to compute the Holevo capacity of a quantum channel. It is a natural but non- 
trivial extension of the existing algorithm for one qubit states. The merit of our 
algorithm is robustness of computation as a global optimization. The algorithm 
by Osawa and Nagaoka [92j can converges to a local optimum and might need 
some iterations of optimization process. Our algorithm overcomes that problem. 
Approximating a continuous object by a point mesh, our algorithm can compute 
a global optimum although it only yields an approximation and its preciseness 
depends on the fineness of the mesh. 

1.4 Organization of this dissertation 

This dissertation is organized as follows. In Chapters [2] and [3], we explain some 
preliminary facts about computational geometry and quantum information theory 
respectively. Our contribution is described in Chapters [H [H and [6l 

In Chapter m we show explicit correspondence between primal and dual quan- 
tum state space. We show some coincidences of Voronoi diagram in the space of 
one-qubit space. The main result of this chapter is divided into two parts: about 
the space of pure states and the whole space including mixed states. This chapter 
is based on the papers [60l |62] . 

The similar problem in a higher level system is described in Chapter O Here, 
we also show the correspondence between primal and dual quantum state space. 
The correspondence is less explicit than one-qubit case but mathematically proven 
as for one-qubit. About the coincidence of Voronoi diagrams, because of its com- 
plicatedness, the space of pure states is only analyzed here. This chapter is based 
on the papers [HI |62] . 

In Chapter [U we propose an algorithm to compute the Holevo capacity of a 
quantum channel and experiment it to show it is really useful. It is also proved 
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that Welzl's algorithm to compute the smallest enclosing ball problem is applicable 
to a quantum state space. The idea of those algorithms was mentioned in |62| . 

The summary of all our contribution is described in Chapter [Tj and we also 
explain a perspective of the future research. 
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Chapter 2 



Computational Geometry 



2.1 Voronoi diagrams and Delaunay triangulations 

First we start with the abstract notation of Voronoi diagrams. 

Definition 2.1 (Voronoi diagram). For a given tuple {X, d, P) where X is a metric 
space, 0? is a distance attached to X, and P = {pi}^i is a set of points of X, the 
Voronoi diagram V is defined as 



and each pi is cahed a site {generator, Voronoi vertex). 

When it is necessary to make a distance associated to the diagram clear, we 
denote V hy Vd- 

In the definition above, the space S and the distance d can be arbitrary. The 
Voronoi diagram most commonly used is a Euclidean Voronoi diagram. It is defined 
as follows. 

Definition 2.2 (Euclidean Voronoi diagram). For a set of sites P{pi}iLi C M", 
the Euclidean Voronoi diagram is defined as 




i=l 



(2.1) 




(2.2) 




i=l 



(2.3) 




(2.4) 
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Just by saying "Voronoi diagram," we usually mean the Euclidean Voronoi 
diagram. Intuitively a Voronoi diagram is a diagram of dominance of sites in the 
terms of the distance. In other words, it is the coloring of the space according to 
which site is the nearest. The following example is the application actually wanted 
by the author. 

Example 2.1 (Nearest station problem). There are four subway stations near 
Hongo campus of the University of Tokyo. The campus is so large that if you 
choose a wrong station to access, it takes unnecessarily a long time to walk. If the 
map of campus were colored as a Voronoi diagram regarding stations as sites, you 
can find the nearest station easily by just telling the color of your current location. 

The region dominated by each site is called a Voronoi region ( Voronoi polygon). 
The edge appears in a boundary of a Voronoi region is called a Voronoi edge. As a 
practical implementation of a Voronoi diagram, just knowing the boundary of each 
region is enough. Each Voronoi edge is a part of bisector line (or in general case 
bisector curve) of a certain pair of sites. Thus, computing a Voronoi diagram can 
be just described as deciding which part of bisector curve appears in a diagram. 
Here, note that some bisector curves do not appear at all. 

The dual diagram for a Voronoi diagram is called a Delaunay pretriangulation 
(Delaunay tessellation). The Delaunay pretriangulation and triangulation of a 
Euclidean space is defined as follows. 

Definition 2.3 (Delaunay pretriangulation and triangulation). For a given set of 
sites P = {pi & M."-}fLi, the Delaunay triangulation D is defined as 

D = |e I e is a line segment between pj and pj, where the bisector of {pi,Pj) 

appears as an edge of Voronoi diagram of P| (2.5) 

If Delaunay pretriangulation is a triangulation, i.e. if for any x E X, there is a 
triple of line segments (si, 52, S3) {si € D) which formulate a triangle and it has x 
as a inner point, then D is called Delaunay triangulation. 

If D is not a triangulation, the triangulation made by joining the vertices of 
non-triangle polygon of D is called Delaunay triangulation 
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In other words, Delaunay pretriangulation is a diagram obtained by connecting 
by a line each pair of sites whose Voronoi region is adjacent. Although definition 
above is only for a Euclidean space, but it can be extended for a general space; for 
a general space, just take a segment of a geodesic curve instead of a line segment. 

2.2 Computation of Voronoi diagrams 

There are some algorithms known for construction of Voronoi diagrams. Here, 
we introduce the incremental method, which is an intuitive algorithm for two- 
dimensinal Euclidean Voronoi diagram. We only explain the outline of the algo- 
rithm to show the main idea. For the detail of the algorithm, refer to [87J. Note 
that the internal data structure to represent the Voronoi diagram is not trivial at 
all, but we skip it. 

The incremental method to construct a Voronoi diagram is described as follows. 
Suppose that a set of point P = {pi, . . . ,pn} is given, we construct a Voronoi 
diagram Vi of the point set {pi, . . . ,pi} step by step and finally obtain the required 
diagram V = V^. Suppose that we are going to add a point p as in Fig. 12.21 The 
rough sketch of the process to add this point is described as follows: 

1. Find a region which the new point Pnew belongs to 

2. Draw a bisector between the new point p and the site which dominate the 
region found 

3. Find neighboring region to the current region 

4. Go back to [2] until it comes back to the original region 

In Fig. 12.21 the process above is explained as follow. First the region by the 
P3 is found and draw a bisector between pi and Pncw The intersection of the 
bisector and the p^'s region appears as a new Voronoi edge ei. Then find a point 
the bisector and the edge of the region meet, and go into the neighboring region, 
the ps's region. Now 62 is drawn as the intersection of the bisector of Pnewi P5 and 
the ps's region. The same process is iterated until it get to the original region, and 
63, 64 and 65 are drawn. 
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Figure 2.1: An explanation of the incremental algorithm to draw a Voronoi dia- 
gram: This shows the situation when Pnew is added to an existing digram formed 
by the sites pi , . . . , ps 

In Fig. 12.21 mathematical rigid computation will guarantee that the start point 
of ei and the end point of 65 is exactly the same. However, in a computer, coor- 
dinates are usually expressed in floating point number, and have some numerical 
error. It means that in Fig. 12.21 the start point of ei and the end point of 65 might 
be slightly different. In usual case as Fig 12.21 it is iiot a problem at all because the 
algorithm only have to maintain the neighboring structure of regions. 

For a general dimensional space, a Voronoi diagram can be computed via a 
lower envelope. Actually, a Voronoi diagram in d-dimensional space can obtained as 
follows. Consider a paraboloid in (i+ 1-dimensional space expressed by x^+i = xf + 
• • and tangent planes at the points which are obtained by lift-up of sites. Then, 
the lower envelope of the tangent planes is a Voronoi diagram (Fig. 12. 2p . Here, a 
lower envelope means the lowest part of a given set of surfaces, and its computation 
is that of convex hull of a polytope. Thus, the computational complexity of a 
Voronoi diagram in d-dimensional space is the same as that of a convex hull in 
d + 1-dimensional space. 

In a (i-dimensional space, the complexity for computation of a convex hull has 



22 



Figure 2.2: An exapmle of a Euclidean Voronoi diagram determined by a lower 
envelope 



been proven to be 0(n log n+riL'^/^J) by several different algorithms |96 tll03ll30lll5j . 
Consequently, the complexity of d-dimensional Voronoi diagram is 0(n log n + 
7^L(f^+i)/2j Note that although this is a polynomial for a fixed dimension d, it is 
exponential about d. 

For a general distance function, a Voronoi diagram can also be considered as 
a projection of a lower envelope of some potential function. Halperin and Sharir 
|43[ m] showed that in three dimensional space, a lower envelope of algebraic 
surfaces can be computed in 0(n^+'^)-time for any small e. Sharir [107] extend its 
result to a d-dimensional space and showed it can be computed in ©(n"^"^"^^). This 
means a Voronoi diagram in d dimensional space whose edges are expressed by 
algebraic equation can be computed in 0{n'^~^^). Another non-Euclidean specific 
distance is analyzed by Icking and Ma [55] . 

2.3 Smallest enclosing ball problem 

The smallest enclosing ball problem is described as follows: for a given set of points 
P, compute the smallest ball which includes all the points of P. The first practical 
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algorithm for the smahest enclosing ball problem is given by Welzl [121J. The 
extension of it to the smallest enclosing ball of balls (SEBB) problem is shown by 
Fischer and Gartner [34], and its implementation is freely available as a part of 
CGAL package [19]. However, SEBB solver in CGAL only works for Euclidean 
distance, and it does not fit our objective. 

The following is Welzl's algorithm to compute the smallest enclosing ball. 

Algorithm 2.1. (Welzl [TU] ) 

procedure minball(P : set of points) i> Compute the smallest enclosing ball 

b_minball(P, 0) 
end procedure 

procedure B_MiNBALL(P,i?) t> Compute the ball which includes P and has R 
in its boundary 

if P = 0ori? = (i+1 (where d is a dimension of the space) then 

return the ball which has R on its boundary 
else 

Choose p £ P 

B ^ b_minball(P - {p} , R) 
ii p ^ B then 

B ^ b_minballP — {p} , R U {p} 
end if 
return B 
end if 
end procedure 

The function b_minball is the main part of this algorithm. b_minball(P, R)) 
computes the smallest ball that includes P under the constraint that all the points 
or R must be on its boundary. This algorithm is based on the idea that more 
constraints on the boundary make the computation easier. So, the chosen point 
p is not included in the current optimal ball, it tries to a new bigger optimal ball 
which has p on its boundary. It works because of the following lemma. 

Lemma 2.1. For point sets P and R, let SEB(P, P) be the smallest ball which 
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includes P and has R on its boundary. Then, for any p & P 



SEB(P, R) = SEB(P -{p},RU {p}) 



(2.6) 



Proof. See [T2T] . 



□ 



Another important point of this algorithm is the first "if" part. Note that if 
^R = d+1, the optimal ball is uniquely determined. If the function is called when 
P 7^ and = d + 1, it simply returns the ball which has R on its boundary. It 
returns a wrong answer if some point of P is not included in the ball constructed 
by R. In that case, to keep the consistency of the specification of this function, it 
is correct to return "undefined" because there is no ball which includes P and has 
R on its boundary. However, our main objective is implement MINBALL and this 
exceptional behavior of b_minball makes it easier. 

Suppose that "or \R\ = d + 1" of the first "if" condition of Algorithm 12.11 
is omitted, and the second "if" condition is replaced by "i? is defined and p 
i?." Denote this different version of function as b_minball'. b_minball' has a 
consistent specification as itself, but returns the same value if called by minball. If 
B_MINBALL is called from minball with #i? = d+l and P = {p} , p ^ R, then the 
condition for second "if" becomes "true" and B_MINBALL is called with #RU {p}. 
^R is already maximum and so, the returned value becomes "undefined." If once 
a returned value is "undefined" in the depth of the call of b_minball, the returned 
value of the top level also becomes "undefined." It is easily checked by induction. 
However, the value of minball is certainly defined, so it is a contradiction. 

By the observation above, we can say it is the same thing whether it checks if 
= d+l in the first "if" or it checks if D is defined in the second "if." However, 
for a practical performance, the original algorithm is better although it doesn't 
affect the order of computational complexity. 

Welzl showed this algorithm ends in an expected 0{n) time if the dimension is 
fixed, where n is the number of the given points (i.e. n = i^P). Its effectiveness 
is shortly based on the low probability of the recomputation of an optimum; the 
probability for the two if-conditions to be true is very low. For the detail of the 
probabilistic analysis, see |121j . 
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Although Welzl's original algorithm is only for the Euclidean distance, Nielsen 
|80j showed this algorithm is also applicable for Bregman divergence, which is a 
pseudo-distance used in classical information theory. Whether it is also applicable 
for the quantum divergence is important for numerical computation of Holevo 
capacity, and is proved to be true as a mostly straightforward corollary of Nielsen's 
result. The detail of its proof is given in Section [6T2l 
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Chapter 3 



Quantum Infomation Theory 



3.1 Quantum states and their parameterization 

A pure state can be expressed as a state vector. A state vector in a c?-level system 
is defined as 

('^i = E«i(^h Ei«iP = i (3-1) 

i=l i=l 

where (•] is Dirac's brakct notation and means mathematically a complex vector in 
C^. The vector {i\ means the i-th element of the orthogonal basis of d dimensional 
complex vector space. Additionally the definition of a state vector is up to a scalar 
multiplication, i.e. {4>\ = ^^(^ (^1 = Si A are equivalent when there exists a 

scalar 7 € such that = 7/3j for all i. 

A density matrix is representation of some probabilistic distribution of states 
of particles. Mathematically it is defined as follows. 

Definition 3.1. A density matrix p is a complex square matrix which satisfies the 
following three conditions: 

a) Hermitian, i.e. p = p* , 

b) The trace is one, 

c) It is positive semi-definite. 
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Moreover, we denote the space of all density matrices of size d x d hy 5(C'^), and 
we called it a d-level system. Especially when d = 2", a d-level system is also called 
an n-qubit system. When d is obvious, we denote 5(C'^) by S. 

A density matrix can express both a pure and mixed states. The state vector {(f)\ 
correspond to |</>) {(j)\ as a density matrix. Here \4>) means a Hermitian conjugate 
of {(pl as in the convention of Dirac's braket notation. A mixed state is a state 
which is not pure. Namely a mixed state corresponds to a state which is mixture 
of multiple states. Actually a density matrix is expressed as: 

d 

P = ^ fflj 1^) (*| di > 0, € M, (3.2) 

i=l 

and the condition for p to be mixed is equivalent to that at least two of Oj are 
non-zero. 

We give another mathematically simple definition for pureness and mixedness 
and also for faithfulness. 

Definition 3.2. A density matrix p is called pure if rank p = 1, mixed if it is not 
pure, i.e. rank p > 1, and faithful if rank p = dim S. 

We also use a notation for subspace of S as follows. 

Definition 3.3. For a given quantum state space S, denote S"^^^^ by 

5P"^ = {p\peS, pis pure} . (3.3) 

^faithful g^j^^ ^nonfaithful g^.^ defined similarly. (Note that especially for one-qubit 

system, ^nonfaithful ^ ^pure^ 

Note that while a quantum state is either pure or mixed, any faithful state is 
mixed. The geometric image of pureness, mixedness, and faithfulness is illustrated 
in Fig. ED 

Especially in two-level system, which is often called one-qubit system, the con- 
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The whole space of quantum states 
(convex object includmg its boundary) 




Figure 3.1: An explanation of pureness, mixedness, and faithfulness 



ditions a), b) and c) in Definition 13.11 are equivalently expressed as 

/ 1 + z X -iy \ 
2 2 

X + iy 1 — z 



x^ +y'^ + z"^ <l, x,y,zG]R. 



(3.4) 



This is called a Bloch ball because it is a ball in the xyz coordinate system. The 
parameterized matrix corresponds to the conditions a) and b), and the inequality 
corresponds to the condition c). 

There have been some attempts to extend this Bloch ball expression to a higher 
level system. A matrix which satisfies only first two conditions, Hermitianness and 
unity of its trace, is expressed as: 
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(3.5) 



Actually, any matrix which is Hermitian and whose trace is one is expressed this 
way with some adequate {6}. This condition doesn't contain a consideration for 
a semi-positivity. To add the condition for a semi-positivity, it is not simple as 
in one-qubit case, and we have to consider complicated inequalities [THIEU]. Note 
that this is not the only way to parameterize all the density matrices, but it is 
reasonably natural way because it is natural extension of one-qubit case and has a 
special symmetry. 

Additionally our interest is a pure state. A pure state is expressed by a density 
matrix whose rank is one. A density matrix which is not pure is called a mixed 
state. A pure state has a special meaning in quantum information theory and also 
has a geometrically special meaning because it is on the boundary of the convex 
object. In one-qubit case, the condition for p to be pure is 



This is a surface of a Bloch ball. On the other hand, in general case, the condition 
for pureness is again expressed by complicated inequalities. 

3.2 Distances 

Before explaining about distances for quantum states, we prepare a mathematical 
notation. 



+ ?/2 + 2^ = 1. 



(3.6) 
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Definition 3.4. Suppose that the matrix p is diagonahzed as 



p = U 



A2 



u* 



V 

with a unitary matrix U . For a given function / 
p is defined as 

/^/(Ai) \ 
/(A2) 



f{p) = u 



V 



u* 



Especially we define 
/VA7 

Vp = u 



\ 



XdJ 



and 



/ log Ai 



log p = U 



log A2 



U*. 



(3.7) 



I, taking / of the matrix 



(3.8) 



(3.9) 



(3.f0) 



V log Ad/ 

Now we define two distances. 

Definition 3.5 (See [l5]). For two pure states p and a, the Fubini-Study distance 
dpsip^o') is defined as 



cos dps {p, 0-) = V Tr {pa) , < dps (p, ^) ^ 2 



(3.11) 



Definition 3.6. For two arbitrary quantum states (i.e. mixed or pure states) p 
and (7, the Bures distance dsip, cr) [16] is defined as 



dsip, 0-) = V 1 ~ '^^ V V^pV^- 



(3.12) 
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Especially if p and a are pure states and expressed as \tf) and jV') respectively, 
the Bures distance is as follows: 

dB(p, a) = ^1 - TV {pa) (3.13) 
= ^l-\{p\a)\ (3.14) 

This means the Bures distance and the Fubini-Study distance are fundamentally 
the same thing for pure states. Moreover, since a space of pure states is defined as 
a unit ball divided by a multiplication, those distances are both natural and have 
a Euclidean-like property. 

3.3 Divergence 

In a classical context, Kullback-Leibler divergence means, in a sense, a "distance" 
from a probabilistic distributions to another. When two probabilistic distributions 
Pi and qi are given, the Kullback-Leibler divergence is defined as: 

DkMq) = Y.Pi\og^ (3.15) 
i ^» 

The quantum divergence is the quantum version of Kullback-Leibler divergence. 
Just like Kullback-Leibler divergence has an important role in classical information 
theory, the quantum divergence is essential in quantum information theory. It is 
defined by a similar formula. 

Definition 3.7. Suppose that two quantum states p and a are given and a is 
faithful. The quantum divergence is defined as 

L>((7||p) =TVa(loga-logp). (3.16) 

Note that though this has some distance-like properties, it is not commutative, 
i.e. D{a\\p) / D{p\\a). Also note that p does not necessarily need to be faithful 
because log can be naturally defined as 0. 

3.4 Quantum channel and its capacity 

A quantum channel is the linear transform that maps quantum states to quantum 
states. In other words, a linear transform F : M(C; cZ) — M(C; d) is a quantum 
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channel if r(5(C'^)) C S{C'^). To preserve the condition for density matrix, there 
is a natural restriction for a quantum channel. A quantum channel T satisfies 
following condition: 

1. it must be trace-preserving, i.e. Tr T(p) = Tr p, and 

2. it must be completely positive, i.e. For any identity map /, the map F®/ maps 
a semi-positive Hermitian matrix into a semi-positive Hermitian matrix. 

Such a map can be shortly denoted by a "TPCP map." In other words, the 
condition for a linear transform to be a quantum channel is to be a TPCP map. 

The Holevo capacity is considered as a classical information capacity of a given 
quantum channel under the consumption that the input state is not entangled and 
the output particle is properly measured. 

Definition 3.8 (Holevo capacity [M!)- The Holevo capacity of a given channel F 
is defined as follows: 



C(F) = max S(^p,F(/>,)) + ^p,5(F(p,)), (3.17) 

1=1 1=1 

where (•) means von Neumann entropy, i.e. S{p) = —plogp. 

Another formulation of the Holevo capacity is given by the following theorem. 
Theorem 3.1 (Ohya, Petz, Watanabe [86j). 

C7(F)= min max D(F(p)||F(a)). (3.18) 

o-e5(C<*) pe5(Cd) 

This theorem means that the Holevo capacity is equal to the radius of the 
smallest enclosing ball with respect to the quantum divergence. From now on, we 
mainly use this smallest-enclosing-ball formulation. 



3.5 Calculation of Holevo capacity 

Our first motivation to investigate a Voronoi diagram in quantum states is the 
numerical calculation of the Holevo capacity for one-qubit quantum states [93] . We 
explain its method in this section. In order to calculate the Holevo capacity, some 
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points are plotted in the source of channel, and it is assumed that just thinking 
of the images of plotted points is enough for approximation. Actually, the Holevo 
capacity is reasonably approximated taking the smallest enclosing ball of the images 
of the points. More precisely, the procedure for the approximation is the following: 

1. Plot equally distributed points on the Bloch ball which is the source of the 
channel in problem. 

2. Map all the plotted points by the channel. 

3. Compute the smallest enclosing ball of the image with respect to the diver- 
gence. Its radius is the Holevo capacity. 

In this procedure, Step 3 uses a farthest Voronoi diagram. That is the essential 
part to make this algorithm effective because Voronoi diagram is the known fastest 
tool to seek a center of a smallest enclosing ball of points. 

However, when you think about the effectiveness of this algorithm, there might 
arise a question about its reasonableness. Since the Euclidean distance and the 
divergence are completely different, Euclideanly uniform points are not necessarily 
uniform with respect to the divergence. We gave partial answer to that problem by 
Theorem 14.51 At least, on the surface of the Bloch ball, the coincidence of Voronoi 
diagrams implies that the uniformness of points with respect to Euclidean distance 
is equivalent to the uniformness with respect to the divergence. 

3.6 Entanglement and additivity problem 

The additivity of Quantum channel is simply stated as follows. 

Conjecture 3.1. For any two channels Fi : 5i ^ Si and F2 : ^2 — > S2, it is 
conjectured that the following equation holds. 



The right hand side of Equation 13.191 means the total capacity when two chan- 
nels Fi and F2 are used separatedly. The difference between Si ^2 and Si x ^2 




(3.19) 
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is entangled states; more precisely, for an entangled state is an element of the set 
E defined as: 

^ = 5i®cS2-{pi®/J2 I Pi G5i,P2 g52}. (3.20) 

So, the additivity conjectme of the Holevo capacity is simply stated as "Do entan- 
gled states not contribute to the capacity of a product channel?" or "Are entangled 
states powerless in the terms of Holevo capacity?" 

The conjecture above is stated in the terms of capacity, but there are some 
direct measures of how much a quantum state is entangled, and similar properties 
for them are conjectured. Now we define some measures. 

Definition 3.9 (Entanglement of formation). For a state p in a bipartite system 
Sa <Sb, the entanglement of formation Ep is defined as 

Ef{p)= ^ min Y.P^^^^ B\v^) {v^\) (3.21) 

where the minimization is over all possible expressions such that p = ^^Pi \ vi) {vi\ 
and J2iPi = 1- 

Definition 3.10 (Minimum output entropy). For a given channel F : 5 — > 5, the 
minimal output entropy 5 of F is defined as 

S{T) = mm S{T{p)). (3.22) 

pes 

The very important fact we have to mention here is that globally, some conjec- 
tured properties about the measures defined above are equivalent. It is stated as 
follows. 

Theorem 3.2 (Shor |109j ). The following four propositions are equivalent. 

1. The additivity of the minimum entropy output of a quantum channel. Suppose 
that two channels Fi and F2 are given. Then 

5(Fi®F2) = 5(Fi) + 5(F2). (3.23) 

2. The additivity of the Holevo capacity ( See Conjecture \3. ij) 
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3. The additivity of the entanglement of formation. Suppose that two states 

Pi G Sai (8i Sbi and p2 G Sa2 <S)Sb2- Then 

Ef{pi ^ P2) = Ef{pi) + Ef{p2), (3.24) 

where Ep is calculated over the bipartite A-B partition. 

4- The strong superadditivity of the entanglement of formation. Suppose a state 
p G Sai 'Si Sa2 S Sbi S Sb2 is given. Then 

Ef{p) > EpiTri p) + Ef{T,2 p), (3.25) 

where Ep 's are over the bipartite A-B system, and Tij means a trace out for 
the space Sai S Sbi ■ 

This means that some conjectures about the measures of entanglement are 
equivalent. All the conjectures are about the power of entangled states. 
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Chapter 4 



Voronoi diagrams for one-qubit quantum 
states and Its Application 

4.1 Overview 

In this chapter, we investigate the one-qubit quantum state space. In this case, the 
theoretic analysis is much simpler than general case because the whole space is a 
three dimensional ball {Block ball). In spite of its simple structure, the divergence 
defined here is still distorted and far from the intuition which we have for a normal 
"distance." 

In this chapter, we show some theorems concerning the coincidences of Voronoi 
diagrams with respect to some pseudo-distances. Although some of the theorems 
are the just special cases of the theorems for a higher level case which we explain 
later, we introduce them here because historically they were proved earlier and the 
proofs are much simpler than the general cases. 

Generally, the importance of investigation for one-qubit system as a start-up for 
a new stream of a research has been emphasized in quantum information theory. 
The same thing can be said for Voronoi diagrams. More can be known for the 
one-qubit system than a higher-level system because of its simplicity. Although 
the Voronoi sites are restricted to be pure states, some coincidences of diagrams 
in a set of mixed states can be proved for the one-qubit system, while in a higher- 
level system, mixed states are much more complicated and similar thing is still not 
known. 
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Our motivation originally started from the algorithm by Hayashi et al. [l^ 
to compute the Holevo capacity of a quantum channel. Since it uses Voronoi 
diagrams with respect to the quantum divergence, it is considered to be important 
to investigate the structure of the Voronoi diagram with respect to the quantum 
divergence. 

In Hayashi et al.'s algorithm, the key factor to attain the accuracy is an approx- 
imation of a set of pure states by discrete points. Especially, how well-distributed 
points you can obtain determines the accuracy of the computation. The wellness 
here is in terms of the divergence. Our result about coincidences of Voronoi dia- 
grams tells that well distributed points in Euclidean space is also well distributed 
in terms of the divergence. Hence, although the divergence is difficult to deal with, 
we can generate well distributed points appropriate for the computation. 

The refinement of the point set is another considerable application. Suppose 
that you have computed the capacity with some generated approximating points, 
and wish to refine the precision of the result. Then, where should the new points be 
plotted? Voronoi diagrams give the answer to that question. Since some diagrams 
with respect to some pseudo-distances are the same, all you have to think about 
is the Euclidean Voronoi diagram. Just plotting new points on the Voronoi edges 
would be a reasonable refinement. 

Although the algorithm by Hayashi et al. is the only existing application of 
our result so far, it can be considered in a general context. The difficulty of the 
numerical computation in quantum information theory is due to a computation of 
a continuous geometric object. It is essentially different from classical information 
theory. Using Voronoi diagrams can be one of the options to overcome the difficulty. 

4.2 Voronoi diagrams in a quantum state space 

We consider Voronoi diagrams with respect to the divergence. Since the quantum 
divergence -D(-||-) is not symmetric, we can consider two different Voronoi diagrams 
about this measure. 

Definition 4.1. For a given quantum state space S, Voronoi diagrams for a given 
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set of sites S = {ai} are defined as: 



Vd{^) = {V^^ = {peS\ D{p\\ai) < D{p\\aj) for any j}} 
= = {peS\ D{ai\\p) < D{aj\\p) for any j}} 



(4.1) 



(4.2) 



The definition of the coincidence of Voronoi diagrams is natural. Suppose that 
a space X and distances di and d2 on X are given. For a given set {pi G X}, the 
two Voronoi diagrams Va^ and V^j are said to coincide when they are equal as a 
set. This property is equivalent to the coincidence of bisector curves. 

To show that given Voronoi diagrams coincide, it is sufficient to check their 
bisector curves. This fact is stated as follows: 

Theorem 4.1. Suppose that a space X and distances di and d2 on X are given 
and they satisfy di{x,x) = for any x G X. Then, the following two conditions 
are equivalent. 

1. For any set of points P = {pi, . . . Voronoi diagrams V^i(-P) and Vd^iP) 
are equivalent. 

2. For any given pair of points {p\,P2), the bisector curve ofpi,p2 with respect 
to di and d2 are equivalent. 

Proof. We show this by induction. For two sites, the only edge is the bisector, and 
because di{x,x) = 0, the region dominated by a site is the same side as the site. 

Suppose that V^j (P) and V^j (-P) are the same for first n — 1 points of P. If 
we add an n-th point as a new site, newly appearing edges are the bisectors of the 
new site and other site. The region dominated by the new site is the same side as 
the new site. Then, the diagrams for n points are also the same. □ 

4.3 Primal and dual Voronoi diagrams 

In this section, we characterize the Bloch ball by means of computational geometry. 
We define primal and dual space of quantum state. We show one of the two 
divergence Voronoi diagrams is linear in the primal space and the other is linear 
in the dual space. All the theorems shown in this section can be extended to an 
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arbitrary quantum state space. However, more explicit computation is possible for 
one-qubit states, and it helps us to understand the things more deeper. It is a 
common merit for an investigation of one-qubit states. 

Theorem 4.2. For one-qubit states, Vd is linear. 

The following is the essential for the proof of this theorem. 

Lemma 4.1. For a one-qubit mixed state a = a{x,y,z) and a general one-qubit 
state p = p{x, y, z), 



r 1. l-r2 



DipWa) 



1 + r 1 , 1 — f 



^2 



■ log 1 — log log ■ 

2^ 4 2^1-r 2^ 4 

1 , 1 + r , _ . 

-— log r [xx + yy + zz) 

2r I — r 



1 . 1 - 



log ■ 
12^4 



1 + r 1 



2 1 — r 



log — 
2 ^4 



(0,0,0)) 

((5,y,5) = (0,0,0)) 

(4.3) 



where the paranieterizations are given as Formula ^3.4^ and r = x"^ + y"^ + z"^ . 
(Note that the formula for {x,y,z) = (0,0,0) is the limit of the formula for 
ix,y,z) / (0,0,0) as (x,y,z) ^ (0,0,0).; 



Note for notation: By a = a{x,y,z), we mean a is parameterized by {x,y,z) 



as in Formula (j3.4p and implicitly define f by f = ^/x^ + ]p + 2^ 

Proof. The eigenvalues of p are 

1 + r 1 — r 
2 ' 2 

When (x,y) 7^ (0,0), defining a unitary matrix U as 



^7 



V2 



I X — iy r -\- z x — iy r — z\ 



V 



r — z 



r + z 



p is expressed as 

/1 + r 

a = U 




(4.4) 



(4.5) 



(4.6) 
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Then, 



Tr p log p =Th: U 



'l+r 






1 - r 



u* ■ u 



( l + r 



Tr U 



/1 + r^ l + r 
' log- 



V 



2 ° 2 




Tr 



'l + r l + r 

log 

2 ^2 





1 — r , 1 — r 



1 — r , 1 — r 
log' 



= - log ■ 
2^4 

If (x,?/) = (0,0), z = r and 

-1 + 2 

Tr p log /3 = I ^ 




2 ° 2 

l+r l + r 1— r 1 — r 

log \ log 

2 2 2 2 

1 . 1 - r2 



r l + r 

+ o log:; ■ 

2 1 — r 



1 + z 




log 





1 - z 



l + z^ l + z 1-2:, 1-2 

= log 1 log 

2^2 2^2 

l+r l + r 1— r 1 — r 

= log 1 log . 

2 2 2 2 

This means Equation (j4.7p also holds for (x,y) = (0,0). 
Now, we will consider three cases: 

. (x,y)^(0,0) 

• (x,y) = (0,0) and 5/0 

• ix,y,z) = 

For {x,y) = (0,0), denoting similarly as p by 



a = U 



/ l + r 
2 






1 -r 



U* 
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we obtain 



IV (p log a) 



= Tr 



pxU 



A 1 + f 
Mog^ 



V 







log- 



1 — f 



= Tr 



pU 



( X + iy If + z 



V2 



log 



1 + f 



If — z ^ 1 + f \ 

— ^log— — 



X + iy f — z ^ 1 — f 
' log ■ 



= Tr 



/^logi±^ + ^f^logif^ 



If + z, 1 — f 



x—iy 
f 



V 



2±iaiog|+f 

r ° 1— r 



r—z 
f 



logl±i + f±Ilogi^y 



= Tr- 



/{i + z) (y log Y + ¥ log + (--^yf+'y) log i±f *\ 

^ * (l-^)(l^logi±i + f±llogi-r) + (£±!^^logi±|^ 

( " * " means an unnecessary element for the computation below) 

1 s f f + z ^ 1+f f — 1 — f\ (x — i ti)(.f + iy) , 1 + f 

=- 1 + z ^log^— + ^log^— +^ ^log,^ , 

4 V?" ^ ^/ 4r 1 — r 

1,, ,/f — 1 + f f + z, 1 — f\ + iu) (x — iw) , 1 + f 

+ - 1-z ^log^— + ^log^— +^ ^log- z 

4 2 r 2 / 4r 1 — r 



1 1 — f^ 1 1 + f 
7, log — A 1- tt: log z {xx + yy + zz) . 

2 4 2r 1 — r 



(4.10) 



Thus, 



L>(p|H =TV (plogp)-IV (ploga) 



1 1 - 

— log 

2^4 



+ o lOo , 
2 1 — r 



1 + r 1 , 1 — f 



— log ■ 
2^4 



1 , 1 + f , , 

- — log z [xx + yy + zz) 

2r 1 — r 



(4.11) 



and this completes the proof for {x,y) ^ (0,0). 
Ifx = y = 0, 5 = f and 



Tr (plogd) 
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Tr 



log 



1 + f 



log- 





1 — f 



1 + z 1 + f 1 — z 1 — f 
— - — log — 1 — log — - — 



1 , l-r' 
- log 



z 1 + f 

+ - log- r. 

2 1 — r 



(4.12) 



Here, 



1 , 1 - P 1 , 1 + f , ^ 

o log — ^ ^ 77^ log r [xx + yy + zz) 

2 4 2r 1 — r 



1 , 1 - f ' 

— log 

2^4 



1 



1 + f 



H log 

2f 1 — f 



zz 



1 l-f^ 1 1 + f 

This completes the proof for x = y = 0, z ^ 0. 
When X = y = z = 0, (j4.12p becomes 

- log - , 

2 ^4' 

and this completes the proof for x = y = z = 0. 



(4.13) 



(4.14) 
□ 



Proof of Theorem \4-^ Suppose that p, ui and 02 are parameterized by {x,y,z), 
{xi,yi, zi), (x2, 2/2,-^2) respectively. Using Lemma we obtain 



D{p\\ai) - D{p\\a2) 



1 - f 1 , 1 - f 2 
log — log — - — 



, - 1 , i + h - 1 , 1 + ^2 

X XlTT^^Og- — -X2— log- — 

lr\ 1 — ri 2r2 1 — 7'2 

, - 1 , 1 + ri _ 1 1 + f 2 

y\y\Tr^ log — ^ - yi^^ log - — — 

lr\ 1 — ri 2r2 1 — 7'2 

, - 1 , 1 + ri .1 1 + f 2 

2; zi— log - 22— log 

lr\ 1 — ri 2r2 1 — r2 



This is linear in x^y,z. 



(4.15) 
□ 



Theorem 4.3. Define a transformation from p = p{x, y, z) to p by 



1 



1-r' 



p = - log p + - log ^ 



(4.16) 
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and 



9 rr. . 

u = -^Tr plog/9 

dy 
d 

w = ——TV plog/9 
oz 

Then, p can be expressed as 
p{u,v,w) = 



W U — IV 

u + iv —w 



Proof, p can be diagonalized as 



p = U 



■1 + r 




\ 

1 — r 



U* 



with 



U 



>r — z 



( X — iy jr + z x — iy jr — z\ 



r -\- z 



Then, 



d 

'^^ 'dx ^""^ P'^^^P) 

dr d , ^ ^ . 
= d^-¥r^-^'^''^'^ 



Here, 



|:(TVplogp) = |;|TrC/ 



d_ 
dr 



Tr U 



l+r 




2 







2 / 


l+r 


logii^ 


2 



u* ■ u 



log 



i^logi^ 



= - iTr 
5r \ 



'^logi±^ 








5 /l + r, l + r 1 — r, 1 
— — log — 1 — log — 



dr 

1, l + r 

-log- , 

2 1 — r 
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and 



dr X 
dx r 



(4.23) 



Thus, 



X 1 + r 

u = -—log- . 

2r 1 — r 



(4.24) 



In the similar way, we can obtain 



y 'i- + r 

V = -—log- , 

2r 1 — r 

z . 1 + r 

w = -—log- . 

2r 1 — r 



(4.25) 
(4.26) 



On the other hand. 



1 l-r2 
log/?+ - log — / 



-U 



log— P^— 



1 -r' 



log 



■I 





1 — r 

-IT) 



( X + iy jr + z 

72 



x + iy 

-|-y2 V r 



r — z 



r — z \ 



r + z 



r ) 



( x—iy fr+£ x—iy / r—z 



r— 2 
r 



V 

1 

- log • 1 

2^4 

/r + z^ l + r . r — z, 1 — r 



.V^^^ r iOg 2 V iOg 2 



log o 
r 2 



log 



V 



X -I- , l + r 
log- 



1-r 



X — iy ^ l + r 

log- 

r 1 — r 



r — z, l + r r + 2; , 1 — r 
log — \ log 



2~/ 



1 1-^2 

/ z , l+r 

' log 

2r ^ 1 - r 

X + iy , l + r 



X — iy , 1 + r\ 

-— log- \ 

2r 1 — r 



z . l + r 
Yr^^'^T^r I 



W U — IV 

^u + iv —w 



(4.27) 
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and this completes the proof. 



□ 



Theorem 4.4. For a one-qubit state p and a, we use plain notation for p and tilde 
notation for a, i.e. a and a are parameterized by x,y,z and u,v,w respectively. 
Define D by 



(4.28) 













V 




V 




V 


\w) 




[w) 




[wj 



D{p\\a) = i;ip) - ipia) - 

where {■,■) means the inner product of vectors and 
ip{p) = log (TV (expp)) . 



(4.29) 



Then, D{p\\a) = D{a\\p) and the Voronoi diagram with respect to D with sites as 
second argument of D (denoted by V^) is linear. 



Proof. 



exp p = exp 



= expC/ 



-U 



log 



1+r 
2 




9 





log 


' 2 , 


y 1+r 
= 2 




9 
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1 1 - 
U* + -log^.I 



1 , 1 - 



■ I 



U* 



^^/'exp(-logi±r + ilogl^) 



U 







exp(-logi-^: + ilogi^)^ 



U* 



1 + r 



'1 + r 
1-rJ 



U* 



(4.30) 
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Thus, 



il){p) = logTr U 



1 + r 



log Tr 



log 



'1 + r 
l-rj 

\ 



U* 



1 + r 







'1 + r 
1-rJ 



U*U 



1 — r 
1 + r 



+ 



'1 + r 
1 — r 



1 1 -r2 
= — log . 

By taking square of (j4.24p . (j4.25p and (14.26P : and adding them, we obtain 



2 I 2 I 2 



log 



1 + r 
1 — r 



(4.31) 



(4.32) 



Taking ^ of this formula. 



2u 



dr d 
du dr 
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log 
log 
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log 



du 2(1 — r^) 1 — r ' 
and similarly, 
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log 
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1 + r 
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Thus, we obtain 



(4.33) 

(4.34) 
(4.35) 
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l-r2 



du 
dr 
dv 

\ / 

\dw/ 



2r log 



1 + r 
1 -r 



(4.36) 



Using (j4.3ip and (j4.36p . D{p\\d-) can be expanded as 



1 1 - 1 1 - f2 
l<7) = -- log + - log 



f u — u\ 

V — V 
\w — w J 



, 2r log 
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1 — r 
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1 l_r2 I i_f2 

- — log i — log 

2 ^ 4 2 ^ 4 

1 1 + r 

— log • {xx + yy + zz) + — log 



2r ° 1 - r 
1 i_r2 1 i_f2 
- — log j — log 



- log ■ {x' + y' + z') 

Ir 1 — r 



1, l + r, _ . r ^ 1 + f 

= Dia\\p) (4.38) 

Note here we used Theorem 14. II Additionally, Formula ()4.37p is linear for u,v,w, 
and consequently the equation for the boundary 

D{p\\ai) - D{p\\a2) = (4.39) 
is also linear. □ 



4.4 Voronoi diagrams for one-qubit pure states 

The definition for the Voronoi diagrams with respect to the divergence in the space 
of pure states is not obvious because the divergence D(p\\a) is not defined for pure 
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a. Actually, while D{p\ \a) = Tr p(log p — log a) can be defined when an eigenvalue 
of p equals because log can be naturally defined as 0, it is not defined when 
an eigenvalue of cr is 0. Here we show that this Voronoi diagram of mixed states 
can be extended to pure states. We shall prove that even though the divergence 
D{p\\a) can not be defined when o" is a pure state, the Voronoi edges are naturally 
extended to pure states. In other words, we ean define a Voronoi diagram for pure 
states by taking a natural limit of the diagram for mixed states. When we say "a 
Voronoi diagram with respect to divergence for pure states", it means a diagram 
obtained by taking a limit of a diagram for mixed states. 

Note that when we Voronoi sites are pure states, we assume a certain kind of 
a parameterization for the convergence. Actually the Voronoi diagram depends on 
how the sites converge. We assume that when sites are given as /3i(si), . . . , Pn{sn) 
and they are all pure, the diagram is considered as the limit of the diagram with 
the sites rsi, . . . , rs„ in the Euclidean coordinate system, where < r < 1 and the 
limit is taken for r — > 1. This definition might seem to be unnatural because it 
trying to fix the way of uncertain convergence but we believe it is natural because 
it is symmetric. Note that this definition is only possible for the one-qubit system 
because it has special symmetry, and in general, for a higher-level system, no longer 
exists. 

To summarize the facts explained above, we give the following definition. 

Definition 4.2. When a set of Voronoi sites S = {ai} is given, the Voronoi 
diagram with respect to the divergence in the space of pure states is defined as 



yP""-^(S) = Closure(vi'^) n <SP"'^^ 
F^P""'"(S) = Closure(V;^'^) n 5P"'^^ 



(4.40) 
(4.41) 



where 



V^^ = lim{p G 5 



a— »1 



'faithful 



D{ai{a)\\p) < -0(^^(0)11^) for any j} 



V, 



D{p\\ 



o-j) < D{p\\aj) for any j}. 



(4.42) 
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Here Closure(-) means a topological closure and ai{a) is defined by 

/ 1 + azi axi — iayi \ 



ai{a) = 

axi + iayi 1 — cLZi 



(4.43) 





(4.44) 



Theorem 4.5. For given one-qubit pure states as sites, the following four Voronoi 
diagrams are equivalent in the space of pure states: 

1. the Voronoi diagram with respect to the Fuhini- Study distance 

2. the Voronoi diagram with respect to the Bures distance 

3. the Voronoi diagram on the sphere with respect to the ordinary geodetic dis- 
tance 

4- the section of the three-dimensional Euclidean Voronoi diagram with the 
sphere and 

5. the two Voronoi diagram with respect to the divergences, i.e. V^^^ and V^^'^^ . 

In Theorem the coincidence of 1-4 is easy to prove. Actually, for p = 
p{x,y,z) and a = a{x,y,z), 



Tr (pa) 



l-\-zx — iyx + iy 1 + zx + iyx — iy 



2 2 2 2 2 2 

1 + XX + yy + zz 



(4.45) 



2 



Thus 



dFs(p,cr) 



COS (Tr (per)) 




(4.46) 



dB(/c,cr) 



- Tr pa 




(4.47) 
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On the other hand, the Euchdean distance is computed as 
dE{p,cr) = \J{x- xY + {y - yY + {z - zf 



{r'^ + P) -2{xx + yy + zz). (4.48) 
EspeciaUy when p and a are pure, since r = 1 and f = 1, 



dE{p,a) = ^2{l- XX- yy- zz). (4.49) 
Then, the equivalence is proved by simple workout as 

'^Fs(/Ol,cr) = f^Fs(P2,0") 



_i I + xix + yiy + ziz _i 1 + X2X + y2y + Z2Z 
^cos y =cos y 

<^ xix + yiy + ziz = X2X + y2y + Z2Z, (4.50) 
and similarly we can easily show 

dB{pi,CF) = dB{p2,cr) 

<^ xix + yiy + ziz = X2X + y2y + Z2Z, 
dE{pi,cr) = dB{p2,cr) 

<^ xix + yiy + ziz = X2X + y2y + Z2Z. (4.51) 

This this complete the proof of 1-4. 

The rest of Theorem 14.51 is proved using Lemma 14.11 and the following lemma. 

Lemma 4.2. For a mixed state a and an arbitrary p, 

D{pi\\a) = D{p2\\a) 

xix + yiy + ziz = X2X + y2y + Z2Z. (4-52) 

Moreover, under the condition fi = f2, 

D{p\\ai) = D{p\\a2) 

<^ xxi + yyi + zzi = xx2 + yy2 + zz2. (4.53) 
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Proof. Because fi = r2, by Lemma HTH 

Dip\\ai) - D{p\\a2) 
1 + fi 

= - log [(^^1 + + ^^i) ~ + yy2 + ZZ2)] ■ (4.54) 

Thus, 

D{p\\a^) = D{p\\<j2) 

^ xxi + yyi + zzi = xx2 + yy2 + (4.55) 

□ 

Note that the setting for Lemma 14.21 is more general than needed for the proof 
of Theorem l4. 51 because p is not restricted to a pure state. Actually, the coincidence 
of Euclidean Voronoi and divergence- Voronoi also holds in the set of mixed states. 
It is proved in the next section. 

Theorem 14.51 means that all the diagrams are the same as the ordinal Euclidean 
one, which has been researched enough. It tells us the computational complexity 
of the diagrams stated as follows: 

Corollary 4.1. In the space of pure one-qubit states, the following Voronoi dia- 
grams can be computed in 0{nlogn)-time for n sites. 

1. the Voronoi diagram with respect to the Fuhini- Study distance 

2. the Voronoi diagram with respect to the Bures distance 

3. the Voronoi diagram on the sphere with respect to the ordinary geodetic dis- 
tance 

4. the section of the three-dimensional Euclidean Voronoi diagram with the 
sphere and 

5. the Voronoi diagram with respect to the divergences, i.e. Vd and Vd*. 

Proof. The geodesic Voronoi diagram on a sphere is computed in 0(n log n)-time 
[TOO]. Then, apply Theorem □ 
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4.5 Voronoi diagrams for one-qubit mixed states 

In Theorem 14.51 we showed some Voronoi diagrams are the same in the set of pure 
states. Some part of the theorem can be extended to mixed states. The coincidence 
of Voronoi diagrams in a set of mixed states is stated as fohows. 

Theorem 4.6. For given one-qubit pure states as sites, the following four Voronoi 
diagrams are equivalent in the whole space S: 

1. the Voronoi diagram with respect to the Bures distance 

2. the Voronoi diagram with respect to the Euclidean distance 

3. the Voronoi diagram with respect to the divergences, i. e. Vd and Vd* ■ 

(Note that here, the sites are restricted to he pure states while the diagram is con- 
sidered in the whole space.) 

Proof. Suppose that aj (j = 1,2) are given as sites. Because they are pure, fi = 
f2 = 1. Under that condition, by Lemma 14.21 



D{ai\\p) = Dia2\\p) ^ Dip\\ai) = Dip\\a2) 

xxi + yyi + zzi = XX2 + ym + ZZ2. (4.56) 

On the other hand, 

d^{p,ai) = dE{p,cri) 



<^ y + 1 - 2{xxi -\- yyi + zzi) = r^ + 1 - 2{xx2 + ym + ZZ2) 
xxi + yiji + zzi = XX2 + ym + zz2. (4-57) 

Thus, the equivalence of the divergence Voronoi diagrams and the Euclidean 
Voronoi diagram is proved. 

Now suppose that the sites aj. {j = 1,2) are expressed as 

(7j = {iPjl (4.58) 

with 

i^j (4.59) 
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which means 

Si^i s,t,\ ^^^^^^ 



Then, since ^Joj = aj the Bures distance for an arbitrary quantum state p is 
calculated as 



dB{crj,p) 
= 1 — Tr y/ojpoj 
= 1 - Tr (ipjlpli^j) 

= i-TVy'|V^,)[(V^,-|p|V^,)] (^,1 



/I + z _ X — iy _ X + iy - 1 — z -A,,,,,, 
= 1 - Tr W iy-^SjSj + -^—tjSj^—Sjtj + -^tjtjj (V'il 

1 + z _ X — iy _ X + iy - 1 — z - , , r.-,\ 
= 1 - W -^SjSj + -^tjSj^—Sjtj + -^tjtj. (4.61) 

Here we supposed p is parameterized as ()3.4p . If fij's are also parameterized in the 
same way as 

2 2 
2 2 ' 

Formula (j4.6ip can be more simplified as 



1 ^j+Wj i-^j ' (4.62) 



dB{(^j,p) = 1 - \j]^{xjX + yjy + z^z) (4.63) 

Thus, 

dB(o-i,p) - dB(o-2,/o) = 

<^=^> (xix + yiy + ziz) - {x2X + y2y + Z2z) = 0. (4.64) 

This means the Bures diagram is the same as the Euclidean Voronoi diagram in 
the xyz-space. □ 

Similarly as for pure states, computational complexity can be known for the 
diagrams: 
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Corollary 4.2. In the space of general one-qubit states, when n sites are given as 
pure states, each of the following Voronoi diagrams can he computed in 0{n'^)-time 

1. the Voronoi diagram with respect to the Bures distance 

2. the Voronoi diagram with respect to the Euclidean distance 

3. the Voronoi diagram with respect to the divergences, i. e. Vd and Vd* ■ 

Proof. The Euclidean Voronoi diagram can be computed in 0(n^ )-time (See Sec- 
tion E^]). Then, apply Theorem S^l □ 

4.6 Meaning of the result 

The direct application of the fact proved above for the coincidences of Voronoi 
diagrams is the algorithm by Hayashi et al. |l47j and Oto et al. [931 El] to compute 
the Holevo capacity. 

In the source of a channel, they plotted points so that they are vertices of 
the mesh which is obtained by dividing the sphere equally both latitudinally and 
longitudinally. It is a intuitive way and looks reasonably well-distributed, but its 
uniformness of the points is in terms of Euclidean distance. In the algorithm, the 
smallest enclosing ball of the images of the points is computed, and it is in terms 
of the "divergence." The accuracy depends on how well the points are distributed 
with respect to the divergence. However, as we have shown, the Voronoi diagrams 
with respect to Euclidean distance and the quantum divergence are the same, and 
it guarantee the uniformness of the points. 

In the original paper by Hayashi et al. |47J and Oto et al. j93|, I94j. there is an 
implicit assumption that uniformly distributed points in terms of the Euclidean 
distance is also uniform in the world of the quantum divergence. The coincidences 
of Voronoi diagrams we have shown partly fill its gap. Since the coincidences we 
proved is only for the case that all Voronoi site are pure, we have not showed the 
uniformness preserve in general case. However, if the image of a given channel is 
sufficiently large (i.e. the surface of the image ellipsoid is near to the unit sphere), 
the images of the plotted points have a similar property as pure states because the 
smoothness of pseudo-distances. 
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Another expected application of our result is a refinement of points. Since 
the algorithm by Hayashi et al. [U] approximates a geometry by some points, it 
becomes more accurate if the number of plotted points becomes larger. Then, where 
should the additional points be located? The reasonable answer is the Voronoi 
edges, because the Voronoi edge can be regarded as the set of the farthest points 
from the existing plotted points. In a general case, the problem arises here is 
difficulty of non-Euclidean distance, but thanks to the theorems we proved, the 
Voronoi edges with respect to the divergence is the same as those of Euclidean 
Voronoi diagram. You can easily refine the point set by adding points on the 
Voronoi edges, not worrying about the distortion of the quantum divergence. 

From the viewpoint of computational geometry, we showed Voronoi diagrams 
can be a tool to compare some measures defined for the same set. The important 
point is that the coincidence of Voronoi diagrams can be a hint for uniformness of 
a point set in different measures, and we showed an example in three dimensional 
space, which can be visualized and observed intuitively. 

4.7 Summary of this chapter 

We proved the coincidences of some Voronoi diagrams for one-qubit states. Some 
part of the result is just a special case of the general result explained in Chapter [H 
but some are specific for one-qubit system. 

More precisely, we showed that when Voronoi sites are given, the Euclidean 
Voronoi diagram V^^, the Bures Voronoi diagram V^^ and the divergence Voronoi 
diagrams Vdi^d* are the same. Additionally, if it is restricted to the set of 
pure states, the diagrams above, the Fubini-Study Voronoi diagram V^^g, and the 
geodesic Voronoi diagram are all the same. Table 14.11 shows the summary of the 
proved facts. 

The known application of it is the algorithm by Hayashi et al. to compute 
the Holevo capacity. It supports the effectiveness of the algorithm, and we also 
suggested the refinement of the point set by add point in its Voronoi edges. 
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Table 4.1: Coincidences of Voronoi diagrams for one-qubit states and their com- 
putational complexities: Note that the Voronoi sites are given as pure states in all 
the cases 



space 


pseudo-distance 


coincidence 


complexity 




Fubini-Study 










Bures 








pure states 


geodesic 




> coincide 


0(nlog n) 




Euclidean 










divergence 


> 








Bures 


> 






mixed states 


Euclidean 




> coincide 


0(n2) 




divergence 


> 
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Chapter 5 



Voronoi diagrams for 3 or higher level 
quantum states 

5.1 Overview 

Theoretic analysis for Voronoi diagrams for three or higher level quantum state 
space is given in this chapter. Our motivation originally comes from the natural 
interest about whether the extension of the theorems shown in Chapter |4] hold or 
not. 

However, the structure of the quantum state space for three or higher level sys- 
tem is much more complicated for one-qubit quantum system. One of the obvious 
support for its complicatedness is the result by Kimura [66], which showed the 
explicit condition for a complex matrix to be a density matrix and the inequalities 
appeared are too much complicated. Consequently, we could not have done analy- 
sis for three or higher level system about the same condition as in one-qubit case. 
We shall only show the case that Voronoi site are given as pure states though the 
Voronoi regions may be general. 

Using our result, we can convert a problem about a certain pseudo-distance 
into another problem about another pseudo-distance. 

5.2 Primal and dual Voronoi diagrams 

The following is the essential property of Vd. 
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Theorem 5.1. The boundaries of the Voronoi diagram Vd are linear. 

Although, only Vd is proved to be linear, the other Voronoi diag ram can be 
obtained by some transformation from a linear Voronoi diagram. It is stated as a 
following theorem; this is based on a common mathematical framework known as 
Legendre transformation. 

Theorem 5.2. Define transformation from p to p by 



p = \ogp+ ixr (logp)/. 



and 



d_ 



Tr (plogp). 



(5.1) 



(5.2) 



Then, the parameterization of p is given as 

m = 



?3d-4 — Csd-S 



Cbd-S + iibd-7 

where p = p{£) is parameterized as in Formula 113. 5|) . 
Moreover, define D by 

D{p\ \a) = ^{p) - ^{a) - (i - fi, V^(77)) , 

where 



did-i - Ef=i ii Cd^-i - 'iS.d'^-i 

^d-l 



(5.3) 



(5.4) 



V'(p) = log(Tr (expp)) . 



(5.5) 



Then, D{p\\a) = D(a\\p) and the Voronoi diagram with respect to D with sites as 
second argument of D (denoted by Vj^) is linear. 
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This theorem is mostly proved by Oto et al. [93], but it contains a logical gap 
or an ambiguous expression. We give a self-consistent proof as follows. 



Proof of Theorem \5.1\ and Denote (p by 

cp = -Tr (plogp), 
and suppose that p is diagonalized as 

p = XAX*, 



A 



A2 



(5.6) 



(5.7) 



(5.8) 



where X is a unitary matrix. Then, 
dip 



dp d\ogp 

. dp \ ( dlogp 

-Tr (^^logpj-Tr [p^ 



(5. 



Here, 



Tr p 



(9 log p 



d 



Tr [p—{X\ogKX*\ 

dX d]np' A 

Tr (XAX*^logAX* +XAX*X-^ ^* 

dX d]np A dX* ^ 

It ( log AX*XAX* — + XKX*X 7' X* + AX*X log A^X 
oCi dii dii / 



X* +XAX*XlogA 



Tr 



+X 



Ad/ 
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TV 



. , . fdx*x . ^ 

AlogA ( ) +X 



Tr 



AIogA|g)+X 



dp 



d 



0. 



\ / ac, Ai 



Ad/ V 



X* 



dKi 1 , 
95, A^/ 



(5.10) 



Thus, 



e. = TV (|^logp). 



(5.11) 



Now compare the each element of Equation (|5.1ip . For z = l,...,(i — 1, only i-th 

dp 



and d-th diagonal elements of ^ are 1/d and —1/d respectively and the other 



elements are zero, i.e., 



dp i-th 



i-th 



l/d 



\ 



(5.12) 



-l/d) 



So, if we write log p element- wise as 



log/3 = (Cij), 



(5.13) 



then 



— J (Cm Cdd)- 



(5.14) 
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Thus, 



d-l 



d-i 



E^« = ^(Ec«-('^-i)Cdd) 



i=l 



i=l 



-Tr log/) - Cdd, 
a 



and 



d-l 



dii -'^ii = Cii- ^Tr log p. 
i=l 



(5.15) 



(5.16) 



Now think of i = d for example, the only non-zero elements of ^ are the (1,2) 
and (2, l)-th, and 



Tr ( ^ log p ) = Tr 



/ 1/2 

1/2 



LV 



log/9 



2C2I + 2C12 

Re C12, 



(5.17) 



because log/9 is also Hermitian. Here Re means a real part. Similarly for i = d+1, 

( -^/2 \ 
i/2 



Tr ( -S^logp] = Tr 



d+l 



V 



logp 



Im C12, 



(5.18) 



where Im means an imaginary part. Similar observation shows Formula (|5.3p . 
In the rest of the proof, we assume 

i=l 
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and it is shown in the lemma below. Then, 

D{piOMr])) = v^{0-^in)-{^-v,i 

= i^iv) - m (e, i) - iv, v)-{^- V, i 

= i'iv) - 4^(0 -{v,'n- 1) 



(5.20) 



We also show 

V^(^) = 7? (5.21) 

as a lemma, and it completes the proof for D{p\\a) = D{a\\p) 

Now we will show both Vd and are linear. For ai and a2, the bisector is 
given as a set of p which satisfies 

D{p\\ai) - D{p\\a2) = 0. (5.22) 

Since 

DipiOMn)) = ^(6 - ^iv) -i^-V, V^(r/)) , (5.23) 
we obtain 

D{p\\ai) - Dip\\a2) 

= ifiC) - vivi) - - vi,^Avi))) - ifiO - ^{m) - - ??2, V^(r?2))) 

= -'f{'ni) + ^{'n2) + V^(ryi)) - (?72, V^(ry2)) - (^,V<^(ryi) - V ^{m)) ■ 

(5.24) 

This is linear equation in p. Because of the similarity of (j5.4p and (j5.23p . Vjj can 
also be proved to be linear in the same way. □ 

Now we show the following lemma to complete the proof of Theorem 15.11 and 
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Lemma 5.1. The following equations hold. 

i=l 



Proof. 



1=1 



Here, 



2 
d 



^3<i-6 + iiid-b 



V" 



d 

id?-2 + ^?d2-l 



Thus, 



i=l 



1 



d 



TV (log p) 



Tr (plog/9) 
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Moreover, 

= ii. (5.30) 

Hence, 

V^{i)=i (5.31) 

□ 

Lemma 5.2. 

The following is another important property of Vd- 
Theorem 5.3. Consider the surface defined by 

C = i^ipiC))- (5.32) 

Then, the Voronoi diagram Vd is obtained as a projection of a lower-envelope of 
tangent planes of this surface at the Voronoi sites (Fig. \ 5.1\) . 



Proof. Because of Formula ()5.23p . the divergence D{p\\a) can be considered as a 
minus the value of the tangent surface at a (Fig. 15. 2p . □ 

Note that this is another intuitive proof of the fact that Vd is linear. Actually, 
since the Voronoi diagram is a lower envelope of planes, its boundaries are linear. 



5.3 Euclidean Voronoi Diagram and divergence Voronoi diagram 

We show that the coincidence which happens in one-qubit case never occurs in a 
higher level case. To show it, it is enough to look at some section of the diagrams 
with some subspace. If the diagrams do not coincide in the section, you can say 
they are different. It is stated as follows: 
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Figure 5.1: An example of a Voronoi diagram obtained from a lower-envelope of 
tangent planes 





c = m 































Figure 5.2: An explanation for a geometric meaning of the divergence 
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Theorem 5.4. Suppose that d> 3 and that the space of general quantum states is 
expressed as Equation h3. 5|) . Then, if given Voronoi sites are in a general position, 
the Voronoi diagram with respect to the Euclidean distance and the Voronoi diagram 
with respect to the quantum divergence V^ are different. 

Note that the diagramed considered here is only. This is because Vd is not 
well-defined for pure sites. Actually, define that 



l + 2e 



0-1 



0"2 



1 + e 



2e 



-de , 



-de , 



0-3 



l + 2e 



(5.33) 



Then, 



-de , 



(5.34) 
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is on the bisector of ai and (73, and is not on the bisector of (T2 and 0-3. However, 

/l \ 



Um (Ti = lim (72 

£^0 £^0 







V 



(5.35) 



This shows the Voronoi diagram depends on how Voronoi sites converge. 
Think the section of p with a d + 1 dimensional plane: 



(5.36) 



Then the section is expressed as: 



d 2 
2 d 







V 







d 



(5.37) 



The elements of this matrix are except diagonal, (0,1), and (1,0) elements. This 
matrix is diagonalized with a unitary matrix as: 



X 

Id-2j 



\ 



|3±1 

d 



X* 

/<i-2; 



where 



d2 



, a +6 +2 , r 

= M + 2' 



A2 = 



6+6+2 r 
2d 2' 



(5.38) 

(5.39) 

(5.40) 
(5.41) 
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2 2 ' 



X 



2d 2 2d 



R. 



4 ^ 2(i 2 



+ 



2d 



(5.42) 

(5.43) 
(5.44) 



Now we will figure out the necessary and sufficient condition for the diagonal 
matrix of Equation (j5.38p to be rank 1. For that condition to hold, the following 
three cases can be considered: 

Case 1 (only d-th raw of the matrix is non-zero) 



6 = 6 



6- 



-1, id = id+i = 0. 



Case 2 (only one i-th raw (3 < i < d — 1) is non-zero) 

6 = 6 = -1,6 = 6+1 = 0, 

all of 6 (3 < i < — 1) are —1 except one (let its index to be k) and = d — 3. 



Case 3 (only A2 is non-zero) 



6 + 6 = - 2, 



6-6 
d? 



+ (d + eJ+i) = i,6 = 6 



6-1 = -1. (5.45) 



Note that it is impossible that only Ai is non-zero. In both Case 1 and Case 2, 
the set of points that satisfies the condition is just one point, so our main interest 
is Case 3. The set of points that satisfies this condition is a manifold. Actually, 
Case 3 satisfies 



(rf_2-2gi)^ 

d2 



+ (d+d+i) = i, 



(5.46) 



and this is an ellipsoid. 
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Then we prepare for workout of the divergence. The log of p is expressed as: 

/logAi \ 



logp 



X 





log A2 

log 13 + 1 



log^^ 



log 



X* 





(5.47) 



Thus, we obtain 



Tr cr log p 



d 



log Ai ^ log A2 



R. 



+ 



+ 



$2-gi I r 



^2—^1 r 



2'^ ^2iogAi + ^^- ^logA2 



R4 



R. 



??2 + 1 



log Ai + ^ — log A2 



R, 



R- 



+ 



With some workout, we get 

R- = r(^-^ + -),R- = -r(^-^-- 
2d 2 r \ 2d 2 



(5.48) 



(5.49) 



Using these fact and the assumption r/i + 772 = 6 + 6 = <^ — 2, we get 
Tr cr log p = 



ed + .d+i6±i^!(!i_-¥) (^^-^ 



2r 



1 ^1 1 1 \ \ 
log — + -logAiA2. 
A2 ^ 

(5.50) 



Next we think of a Voronoi diagram with only two regions for simplicity. It is 
enough for our objective. Let a and a be two sites, and suppose that p moves along 
the boundary of the Voronoi regions. Suppose that a and a are parameterized by 
{r]j} and {fjj} respectively in the same way as p. 

We consider what happens if r(0 < r < 1) is fixed and the following holds: 



6 + 6 = f^-2,6 



6-1 



(5.51) 



70 



The condition < r < 1 means that p is semi-positive and not a pure state while 
r = 1 in pure states. In other words, we regard that p is on the same elhpsoid 
obtained by shrinking the elhpsoid expressed by Equation (|5.46p . These settings 
are in order to take a limit of a diagram to get a diagram in the pure states. Taking 
the limit r — > 1, we can get a condition for pure states. This procedure is analogous 
to the method used in Chapter [H 

Now to think of the shape of boundary, we have to solve the equation 



D{a\\p)=D{d\\p), 
and this is equivalent to 

Tr (c7 - ct) log p = 0. 
Using Equation (j5.50p . we obtain 

Tr (cj — fj) log p = 

{Vd-Vd)U + {Vd+i-Vd+i)U+i 



(5.52) 



(5.53) 



1 

2r 



4(r?i-r?i)(ei-^ 



log^. (5.54) 



Here when r = 0, this is zero because A1/A2 = 1. In that case, p can take only one 
point, but we do not have to care about this case because we are going to take the 
limit r — > 1. From now on, we suppose r > and that means A1/A2 7^ 1. 

Hence we get the following equation that holds in the boundary of the Voronoi 
diagram: 

4(m - m) (6 



ivd - fid)id + {m+i - Vd+i)Cd+i + 



d-2 
2 



0. 



(5.55) 



Consequently, taking the limit r ^ 1, we get Equation (j5.55p as the expression of 
the boundary in pure states. 

A careful inspection of Equation (j5.55p tells us a geometric interpretation of 
this boundary. We obtain the following theorem: 

Theorem 5.5. On the ellipsoid of the pure states which appears in the section with 
the {d+ l)-plain defined above, if transferred by a linear transform which maps the 
ellipsoid to a sphere, the Voronoi diagram with respect to the divergence coincides 
with the one with respect to the geodesic distance. 
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Proof. Think of the affine transform defined by 



y 

V ' 1 



id 



(5.56) 



then Equation (j5.55p is expressed as 

x'{x - x) + y{y -y) + z {z - z) = 0, 
while Equation (j5.46p becomes 



x^ + y^ + z^ = 1. 



(5.57) 



(5.58) 



Thus when (x,y,z) and {x,y,z) are fixed, the point {x',y',z') which stand for r] 



runs along the geodesic. 



□ 



Now we work out the Voronoi diagram with respect to Euclidean distance. 
Under the assumption above, the Euclidean distance is expressed as 

d{(T, p) 

= im - iif +im - 6)^ + iVd - idf + iVd+i - id+i? 

= 2(r?i - ei)2 + {Vd - id? + (%+i - id+i)\ (5.59) 
and we get the equation for boundary as 

d{a,p)-d{a,p) = -4(??i-??i)^i-2(r/d-r?d)^rf-2(r/d+i-??rf+i)^d+i + 2(r/f-7/^) 

+ {vl-fd) + ivl+i-fd+i) = 0. (5.60) 

By comparing the coefficients of ^i, ^rf, and id+i, we can tell that the boundaries 
expressed by Equation (j5.55p and (|5.60p are different. To show how different they 
are, we give some examples in the rest of this section. 

Example 5.1. Suppose that (??i,%,%+i) = {d - 1,0,0) and (rfi, ??d+i) = 
(—1,0,0), then the boundary is = for the both diagrams. 

Example 5.2. Suppose that (??i, %+i) = (0,1,0) and (rfi, ??rf+i) = (0,-1,0), 
then the boundary is, for both the divergence and Euclidean distance, expressed 
by id+i = 0. 
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Example 5.3. Consider the Voronoi diagram with the following eight sites: 



where it's mean all the possible combinations. Then the Voronoi diagrams look 
like Fig. l5.3[ This figure is also for d = 5. Obviously they are different. 



Figure 5.3: An example of a diagram which appears as a section of a Voronoi dia- 
gram in three level quantum state space. The left is the diagram by the divergence, 
and the right is by the Euclidean distance. 

5.4 Other Parameterization 

We have shown that the Euclidean Voronoi diagram and the divergence Voronoi 
diagram are different as far as the regular parameterization of a quantum state 
is used. However, the parameterization of the quantum state is not unique. The 
condition that a matrix is Hermitian and its trace is one is expressed in another 
way. Here, we show an example of the parameterization with which the section of 
the diagrams shown in the previous section coincides. Note that it just means the 
sections are the same and it cannot conclude that the diagrams are globally the 
same. 

It is difficult to investigate all the possible parameterization, or equally all the 
possible embedding of quantum state space into the Euclidean space. This section 




(5.61) 
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is intended only to show an example to which the proof of the previous section 
cannot be applied. 

Suppose that a density matrix is parameterized in another way than Equation 
(j33]) as 









^3d-4 


- i^sd-s 


d 


2 






2 

- i^bd-i 


2 


d 






2 










- i^d-^-i 



V 



6d-4 + 



d-3 



d 

(.d'2-2 + i^d-^-i 



2 



(5.62) 



d 



J 



In other word, we will think what happens if .^j(0 < i < d — 2) is replaced by "^Ci- 
We will show that Voronoi diagrams with respect to Euclidean distance and the 
divergence coincide for pure states in the section expressed by (j5.36p 
Under this parameterization, Equation (I5.46[) is expressed as 



d-2\ 



^/2d 

and Equation (j5.55p becomes 

(% - f}d)id + {Vd+i - Vd+i)^d+i + 2(??i - m) (^6 - ^7^) 



(5.63) 



(5.64) 



Assuming that Equation (j5.63p holds, the Euclidean distance is calculated as 



d{p, a) = 2(r/i - Cif + (% - Cdf + (%+i - Cd+i? 



2 Ui 



d-2\ 



d-2\ 



V2d) ~ ~ '72d) 
d-2\' 
V2d) 



+ Tld + Vd+i 



+ 



+ (% - Cd? + (%+i - ^d+if 



d-2\ 



d-2\ 
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Therefore, 

d{p,a) - d{p,d) = 

- 4(r/i - fji) - ^^-^^ - 2(% - - 2(7?d+i - = 0. 

(5.66) 

This is equivalent to Equation (j5.64p . and we have shown the Voronoi diagrams 
are the same. 

Again, we mention this does not conclude that the diagrams are the same 
globally. We conjecture that they are different. More generally, we conjecture that 
even with any parameterization of quantum state space, the Euclidean Voronoi 
diagram and the divergence Voronoi diagrams are different. 

5.5 Bures distance and Fubini-Study Distance 

In this section, we prove the following theorem: 

Theorem 5.6. In a general level quantum system, for pure states, the following 
diagrams are equivalent: 

• diagram with respect to the divergence, i.e. Closure(VD)* H 5^"''^ 

• diagram with respect to Fubini-Study distance 

• diagram with respect to Bures distance 

The equivalence between the Fubini-Study diagram and the Bures diagram is 
obvious because 

dsiPjCr) < dsip,^) 



\/l — Tr pa < \/l — Tr pa 
Tr pa > Tr pa 

cos^^ \/Tr pa < cos^^ ^Tr pa 
dFs(/0,cr) < dFs{p,&) 

(5.67) 
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Hence we will show the coincidence between the diagram by Bures distance and 
the diagram by divergence. 
For e > G M we define 

fl-{d-l)e \ 



Pe 



X 



(5.68) 



where 



X is a unitary complex matrix expressed hy X = (xij). 



(5.69) 



In other words, is parameterized by one real parameter e and cf complex pa- 
rameters Xij. Note that in this parameterization, there might be more than one 
representations for one given density matrix. 

This Pe converges to a pure state when you take the limit e | 0. Moreover p^ 
has rank d if e > 0, and all pure states in space of d level system can be expressed 
by the limit of p^. Actually, taking the limit, we obtain 



lim Pf 

40 



/xilXii XllXu 
XuXll Xy2X\2 

\X\dX\\ XxdXvi 



X\\Xxd\ 
XuXld 

XldXldJ 



fxu\ 

Xl2 

\xidj 



,xn xi2 



xid) , (5.70) 



and any pure state can be written in this way. 

About the divergence, we will calculate the boundary of the Voronoi diagram 
as a set of p^s for a fixed e, and take its limits. Considering the equivalence of the 
equations stated in (I5.67p . we can restate the problem: 



Claim 5.1. For given pure states ai and a^, let 

5'e(<Tl,Cr2) = |pe D{(Jx\\Pe) = D{a2\\Pe)\ 

and 

T{a\,a2) = \p Tr (fjip) = Tr (fT2/9)|. 
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(5.71) 



(5.72) 



Then, for any (J\^ai the following equation holds: 



lim Se{ai,a2) = T{ai,a2) 

e—^0 



Proof. Suppose that aj's are expressed as 

(m m ■■■ Vid 



yi2 



) (^ = 1,2). 



\yid/ 

Defining tlie matrix Y by 
fyii y2i • 



Y = 



yi2 2/22 



\yid y2d 

we can also express 



0\ 


0/ 



ai=Y 



\ 

Now we obtain 











Y*, = y 



Y* 



D{ai\\p) - Dia^Wp) 

= Tr (cJi log CJi - (Ti log p) - Tr ((72 log (72 - (72 log p) 

=Tr ((72 — (7i) log p (because Tr ai log (7i = Tr (72 log (72 = 0) 



=Try 



/log(l-((i-l)e) 



Y*X 



loge 



(5.73) 



(5.74) 



(5.75) 



(5.76) 



X* 



loge/ 
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=Tr 



/log(l-((i-l)e) 



Y*X 



loge 



Denoting X*Y hy Z = (zij), (j5.77p can be expanded further as 



(I5T71) =Tr 



/log(l-(d-l)e) 



Z* 



log 



=Tr 



Z* 



( zn log(l -{d- l)e) zi2 log(l - (d - l)e) 
221 log e Z22 log e 



V Zdi log e 



^^^2 log e 



=Tr 
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zuzii log(l - {d- l)e) + Z21Z21 log e H h ZdiZdi log e ? 

^12^12 log(l - {d- 1)1) + Z22^22 log e H h Zrf2^d2 log £ 



( "?" stands for a non-zero element that does not 
affect the result of the calculation) 

= - [ziizu log(l - {d- l)e) + Z21Z21 log e H h ZdiZdi log e] 

+ [^12^12 log(l - (d - l)e) + Z22Z22 log e H h 2^22^2 log e]. 

(5.79) 

Here note that Zij = for j > 3 because the elements X* are all zero except the 
first two columns. Thus, we get 

D{ai\\p) - D{a2\\p) = 

|zii|^log(l - {d- l)e) + |2;2iplogeH h jz^i^loge 



+ 



|zi2plog(l - {d- l)e) + 1^221^ log eH h \zd2\^loge 



Zll\ 



+ i F12 

<^=^ (l^l2p - 1^11 

Here we used the fact that 



log(l - (d - l)e) - loge +logej 

log(l - {d- l)e) - log e] + log e }> = 



log(l-((i-l)e) -loge 



0. 



UllP + \z2l\'^ H h UdlP 



kl2|^ + ^221^ H h Ud2p 



1, 



(5.80) 



(5.81) 



because X is unitary and vectors (yn, yi2, . • • , yid) and {1/21, 1122, ■ ■ ■ , y2d) have a 
unit length. For Equation (|5.8Up to hold independent of e, the necessary and 
sufficient condition is 

\zii\ = \zi2\. (5.82) 

Writing down the elements by Xjj's and yjj's, we obtain the condition written as 

l^iiyii + X122/12 H h xidyid\ = 1^112/21 + ^122/22 H h xidy2d\- (5.83) 
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Now we consider the condition for Bures (or Fubini-Study) distance. Similarly 
we can extend the formula as follows: 



Tr o-ipo - Tr a2Po 



--TrY 



1 



=Tr 



=Tr 



-1 



Y*X 



X* 



(i ^ 





V oy 



(i \ 



z 



V V 



fznzn zuzi2 ■■■ znzid\ 

Z\2Z\\ ZX2Z\2 ■■■ Zi2Zid 
\zidZll ZldZu ■■■ ZidZldJ 



--ZuZn - Z12Z12- 



Thus, 



Tr (cji - a2)po = \zu\ = {z^l, 

and finally we obtain the same condition as (j5.83p . 



(5.84) 

(5.85) 
□ 



5.6 Expected applications 

Although the Euclidean distance is an exception and the space is restricted to pure 
states, we have shown that the coincidence of the Voronoi diagrams also happens 
in three or higher system. This means some problem about a distance for pure 
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states can be translated into another problem in another distance. In such a sense, 
we can say we have clarified the structure of pure states as a space. 

As is explained in Section [4.61 effectiveness of the algorithm by Hayashi et al. 
to compute the Holevo capacity of one-qubit channel is partially supported by the 
coincidence of Voronoi diagrams. It is natural to think the extended coincidence for 
higher level system might become useful in the extended algorithm. The extended 
algorithm is introduced in Chapter [H but we have not found a concrete application 
in it. We also have to mention that the main part of the theorem about the 
coincidence used to support the effectiveness of the algorithm for one-qubit system 
is the one between Euclidean Voronoi and the divergence Voronoi, and it is proved 
not to happen in a higher level system. 

The distances in a quantum state space are originally considered as a tool to 
distinguish the states by measuring. We expect our result will be used in the 
identification of states. Actually, how well a message is coded in quantum system 
partially depends on how uniformly distributed points we can get in a quantum 
state space. The uniformness of a point set is preserved between multiple pseudo- 
distances if their Voronoi diagrams are the same. 

As a real vector space, a d-level quantum state space has a quadratic dimension 
d? — 1. Prom the viewpoint of computational geometry, our contribution is an 
implication about such a large dimensional space about its actual computation. It 
provides a new category of research fields: distortion measure in a high dimensional 
space. 

5.7 Summary of this chapter 

We have investigated whether the same thing as the coincidence of Voronoi dia- 
grams in the one-qubit system also occurs in three or higher level system. The 
Euclidean Voronoi diagram and the divergence Voronoi diagram are proved to be 
different in three or higher level system. On the other hand, for pure states, the 
coincidence of the Euclidean distance, the Bures distance, and the Fubini-Study 
distance occurs even in three or higher level system. 

Some problems about a pseudo-distance for pure quantum states can be trans- 
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lated into a problem about another pseudo-distance. There is still no concrete 
application of it, but by the analogy of the numerical computation of capacity of 
one-qubit channel, it is likely to become useful in the future. 

Another point is that we found a connection among measures which have ap- 
parently have no relation. Especially, the connection between a distance used for 
quantum database search and a pseudo-distance used in quantum information the- 
ory is meaningful. 



82 



Chapter 6 



Numerical Computation and Experiment 

6.1 Overview 

We propose a new algorithm to compute the Holevo capacity of a quantum channel 
especially for three or higher level system. Our algorithm is based on a global 
optimization and never converges to a local optimum. It is the main merit of our 
algorithm compared to the preceding one. 

Osawa and Nagaoka [92] proposed an algorithm to compute the Holevo capacity, 
which is an extension of the algorithm by Arimoto and Blahut to compute the 
capacity of a classical channel. The main deficit of the algorithm by Osawa and 
Nagaoka is that it only computes a local optimum, although Arimoto-Blahut's 
algorithm is guaranteed to reach the global optimum. It is because the objective 
function in the optimization is concave for a classical capacity while it is not concave 
for a quantum one. 

The algorithm we propose is an extension of the one by Hayashi et al. [47] and 
Oto et al. 193^194]. We follow their idea of using the smallest enclosing ball problem 
after approximating the geometric object by a set of points. However, to solve the 
smallest enclosing ball problem, the method based on a farthest Voronoi diagram 
explained in [931 IMl cannot be directly applied to a higher level system because of 
its complexity. Our main contribution is to show theoretically and practically that 
Welzl's algorithm to solve the smallest enclosing ball problem is also useful for a 
quantum state space. 
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(a) In Osawa-Nagaoka's algorithm, the 
objective function is not concave and it 
might converge to a local optimum 



(b) In our proposed algorithm, the small- 
est enclosing ball problem is computed for 
plotted points and it is the global opti- 



mum 



Figure 6.1: Conceptual explanation of Osawa-Nagaoka's algorithm and our algo- 
rithm 

The main merit of our algorithm is that it can find a global optimum although 
it is only an approximate algorithm. The difference of the ideas between our 
algorithm and Osawa-Nagaoka's one is conceptually explained by Fig. 16. 1[ 

6.2 Smallest enclosing ball problem in a quantum state space 

In this section, we show that Welzl's algorithm is also applicable for the quantum 
state space with the quantum divergence as its pseudo-distance. 

Matousek and al. [75] showed that the smallest enclosing problem belongs to 
the class called linear programming-type (LP-type) problem which is introduced as 
an abstraction of some algorithm for linear programming. The LP-type problem is 
defined as follows [75] . Think about the class of the optimization problems defined 
by pairs {H, w), where H \s a. finite set and w ^ W is a function with a values 
in a ordered set W . The goal is to find the optimal set Bh which satisfies 

w{Bh) = w{H) and for any G C Bh, w{B) < w{Bh)- 

This problem is LP-type when the following axioms are satisfied. 
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Axiom 6.1 (Monotonicity) . For any F and G with F G G <Z H, we have 



w{F) < w{G) (6.1) 

Axiom 6.2. (LocaHty) For any F C G C H with (— oo <)w{F) = w{G), and any 
heH, 

w{F) < w{F U {h}) imphes w{G) < wiG U {h}). (6.2) 

The smallest enclosing ball problem in Euclidean space is obviously belongs to 
this class. For the smallest enclosing ball problem with respect to a given arbitrary 
distance in a given space, whether Welzl's algorithm is also applicable, or in other 
words, has the same complexity order as in Euclidean space is certified by checking 
if the axioms above are satisfied. 

Nielsen and Nock [80] showed Bregman divergence satisfies the axioms above. 
We will show that the quantum divergence also satisfies the axioms in the rest of 
this section, but it is by mostly the same idea of Nielsen and Nock. 

The latter axiom, locality, can be proved by showing the uniqueness of the 
optimal ball. The uniqueness is follows from the linearity of the bisectors (Theo- 
rem [5?T]). 

Note that there are two types of the smallest enclosing ball problem can be 
considered as 

SEBz5(P) = min maxD(p||a), (6.3) 

pg^faithful o-eP 

SEBUV)= mm maxD{a\\p), (6.4) 

pg^faithful aeV 

where V means a given set of quantum states. Our main interest is the second 
type because it is used to compute the Holevo capacity. The locality can be proved 
for both. First it is directly proved for the first type. For the second type, we take 
dual of the problem as in Theorem 15.21 and similarly in the dual space, we can 
determine the smallest enclosing ball uniquely. 

We have shown that for the smallest enclosing ball problem with respect to 
the quantum divergence, we can reach the optimum following the same process 
as Welzl's algorithm. However, note that it is only under the assumption that 
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the most primitive part of the algorithm works. In other words, it can come to 
no conclusion about the complexity. Actually, it needs a complicated non-linear 
optimization to find a ball in terms of the divergence which passes through given 
points. Not only the complexity for it is not known, but there is no guarantee that 
the solution can be found. 

6.3 Algorithm for numerical computation of the Holevo capacity 

For the fixed level d of quantum system, the top level procedure of the algorithm 
we propose is as follows: 

Algorithm 6.1. 

1: procedure MAlN(r : quantum channel) > Compute the Holevo capacity of a 

given channel 
2: P <— DIST_POINTS 
3: r(P) 
4: B QMINBALL((5) 
5: return radius of B 
6: end procedure 

The function DlST_POlNTS computes the reasonably distributed points to some 
extent, QMINBALL computes the smallest enclosing ball with respect to the quantum 
divergence. As an algorithm for dist_points, we employed the following: 

Algorithm 6.2. 

1: function DlST_POlNTS > Compute the reasonably distributed point in the 
quantum state space 

2: ¥'= (¥'1,...,</C2(d-1)) ^ (0, ...,0) 

3: repeat 

4: V' - 1 - E?'~'^ V 

5: $ ^ {ipi +iip2, ifS+i'Pi, ^2d-3 + iV2d-2, V')^ 

6: R = RU{\^) {^\} 

7: (/P NEXT_STATE((^, 1) 

8: until If is null 
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9: return R 
10: end function 

The function next_state called from this is as follows. Here, A is a constant 
given as a parameter. 

1: function next_state((/9,z) 

2: if I > 2(d- 1) then 

3: return null 

4: end if 

5: (fi ^ Lpi + A 

6: ii ip^ > I - {ifi -\ h (fi^i) then 

7: <Pi^O 

8: return next_state((^, i + 1) 

9: else 
10: return 
11: end if 
12: end function 

This algorithm is equivalent to run through all the possible tuples of inte- 
gers (ni, . . . ,n2(d-i)) which satisfies (riiA)^ + • • • + (n2(rf_.i)A)^ < 1. The similar 
mechanism as in a carry-up in the computation of sum of two integer is used in 
NEXT_STATE. It tries to add A to the i-th value (pi, but if the resulting value does 
not satisfy {ipi, . . . , f2{d-i)) < Ij (/"i is set to and it tries to carry up to the next 
value 

Since |<I>) (<I>| has the same value up to the multiplication of a complex number 
z to ^ (of course, under the condition that \z\ = 1 because \z^) must be a pure 
quantum state), d-th value of vector $ can be restricted to M. So, lines H] and O of 
DiST_POiNTS do not lose generality. 

We are not convinced that this is the best algorithm for dist_points, but we 
take a trade-off between the simpleness of the algorithm and the goodness of its 
distribution. The algorithm is easy to understand, and looks reasonable to generate 
uniformly distributed points. Nevertheless, there is no theoretical guarantee that 
it really generates uniformly distributed points. 
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The main part of our algorithm is the procedure to solve the smallest enclosing 
ball problem. It is described as follows: 

Algorithm 6.3. 



1: procedure qminball(P : set of points) 
2: B_QMINBALL(P, 0) 

3: end procedure 

4: procedure B_QMINBALL(P,i?) 
5: if P = or i? = then 
6: return b_qmb(P) 

7: else 

8: Choose pe P 

9: S B_QMINBALL(P — {p} , R) 

10: if p P then 

11: B b_qminballP — {p} , R U {p} 

12: end if 

13: return B 

14: end if 

15: end procedure 



The hardest part is the algorithm for B_QMB. B_QMB computes the sphere which 
passes through all the given points. To implement B.QMB, the non-linear (convex) 
optimization is necessary, and we have not found the algorithm to compute it 
certainly. However, we believe that for a small d, we can use some general-purpose 
optimization libraries. We show in the next section that it is a really practical 
option. 

Although there is no 100% reliable way to compute b_qmb, the characteristic of 
the algorithm that it runs through points globally is a merit to cover that problem. 
To approximate the final result, it is not necessary to compute the exact value for 
each call of b_qmb, but it is enough with mostly correct value for mostly all the 
call of B_QMB. Since it runs through all the plotted points globally, it is enough 
if the correct value is computed in one of the neighbor of the real optimal point. 
Of course, here we assume the error is one-sided, i.e. when it returns a wrong ball. 
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the ball is always bigger than the real smallest ball. 

The computational error must be considered after divided into two phases. 
About the smallest enclosing ball problem of points, the error is one-sided. In that 
phase, the error is because of the optimization process, and if the optimization 
does not converges to the real optimum, the computed value is always bigger then 
the real optimum. However, our final objective to compute the Holevo capacity 
and it is equivalent to solve the smallest enclosing ball of a continuous geometric 
object. Since a geometric object is approximated by plotting points, one-sidedness 
of the error is not guaranteed at all. 

Of course, we can say we can obtain a better result by plotting more points. 
Then, how many points are needed to achieve a given upper bound of the error? It 
is left to be open. It is essential to make the proposed algorithm really practical. 

6.4 Numerical experiment 

Now, we show the algorithm we proposed is really practical. We will check whether 
the samples computed by Osawa and Nagaoka [92j is correct. 

The most primitive part to compute the sphere with respect to the di- 
vergence which passes through given points (b_qmb in Section 16. 3p is im- 
plemented using GNU Scientific Library (GSL) [42j. We used the function 
gsl_inultimin_f unctionJ df in GSL to compute the minimum value of a given 
function with its domain in R". The optimization problem we have to solve is 
formulated as follows: 

Input: n < d? and ai, 
Minimize D{ai\\p) 
Subject to D{ai\\p) = 

However, GSL has no non-linear optimization function which process the con- 
strains like this. To fit to the specification of gsl_multiminjEunctionjEdf of GSL 
we rewrite the problem with a sufficiently large number A as follows: 



...,an 

■■■ = D{an\\p) 
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Input: n < (P and cii, . . . , ct^, 



Minimize D{ai\\p) 

iD{ai\\p) - DiaMf + ■■■ + - Z?KII/>)) 

(6.5) 

Note that the only exception is the case n = cP. In that case, the computation 
of the center of the bah does not include optimization. The center is uniquely 
determined by equations. Actually, the equation 



DiaiWp) = Dia2\\p) 



DM\p) 



(6.6) 



includes — 1 equal marks and (f — 1 unknown variables, and it has a unique 
solution in general. 

We employed a restricted version of the algorithm which ignore the ball with d? 
points on its boundary. In other words, we solved the problem under the assump- 
tion that the optimal ball is determined by less than points. This restriction 
is because there is no easy way to compute the solution of Equation (16. 6p . Al- 
though this restriction seems unnatural, we expect that it makes no difference for 
the result, or at worst, the difference is small because of the characteristic of the 
problem. 

We compute the capacity for the channel Ts in [92]. Ts is given by 



T{p) = VipV,* + V2pV^ + VspV^* 

where Vi , V2 , V3 are given as 

^0.2 0.3 0.4\ 
Vi= 0.5i , 
^O.li OAi 0.5iy 

^0.1-0.3i 
V2= -0.3i 

y0.3-0.3i 0.2-FO.li 



(6.7) 



\ 
0.1 -0.2« 

/ 



I-ViVi - V^V2. 



(6.8) 

(6.9) 
(6.10) 
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Result by 

Osfiwa and Nagaoka 



0.5 0.25 0.2 0.1 



A 



Figure 6.2: Result of numerical computation for various A 



The step value (A in Algorithm 16. 2p is set to 0.1, and the number of plotted 
points are 49486. This value is determined so that the computation ends in a 
reasonable time, and returns a reasonable result. 

The result of our experiment was 0.672 • • • , while the optimal value written in 
|92j is 0.677 • • • . The arguments to achieve the optimum is shown in Table 16.11 
where the optimal value for p is expressed by 



p* = ^Wiai, 

i 

and the approximated capacity is given by 
C = D{ai\\p*) = D{a2\\p*) = ■ ■ ■ ■ 



(6.11) 



(6.12) 



The machine used in experiment is Intel Xeon 2.3GHz (64bit) with Linux x86_64 
installed. The time for computation was 6 hours and 55 minutes. 

Fig. 16.21 shows the result of the computation for various A. Since we can 
approach to the real solution with a smaller A, it seems our result is approaching 
to the result by Osawa and Nagaoka. However, we cannot conclude because we 
have not estimated the computational error and the error for our algorithm is not 
one-sided. 
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Table 6.1: Result of numerical computation 

Result: 0.6729054 

Time for computation: 6h55m 

Variables to achieve the optimum: 

Wi 



0.082689 



0.15770 



0.13717 



0.17875 



0.21126 



0.23243 



/ 0.65000 + 0.00000i 
-0.37000 + 0.29000i 
\ 0.08000 -O.OlOOOi 
/0.16000 + 0.00000i 
0.28000 - 0.04000i 
^0.23324 + O.OOOOOi 
/0.17000 + 0.00000i 
0.29000 + O.OSOOOi 
(^0.22978 + 0.05745i 
/ 0.01000 + O.OOOOOi 
0.06000 - O.OOOOOi 
\-0.07937+0.00000i 
/ 0.05000 + 0.00000i 
0.10000 - 0.05000i 
y-0.16733 + 0.08367i 
/ 0.73000 + 0.00000i 
-0.36000 + 0.23000i 
»^ 0.11314 + 0.04243i 



-0.37000- 0.29000i 
0.34000 + 0.00000i 

-0.05000 - 0.03000i 
0.28000 + 0.04000i 
0.50000 + O.OOOOOi 
0.40817 + 0.05831i 
0.29000 - 0.03000i 
0.50000 + O.OOOOOi 
0.40212 + 0.05745i 
0.06000 + O.OOOOOi 
0.36000 + O.OOOOOi 

-0.47624 +0.00000i 
0.10000 + 0.05000i 
0.25000 + 0.00000i 

-0.41833 + 0.00000i 

-0.36000- 0.23000i 
0.25000 + 0.00000i 

-0.04243 - 0.05657i 



0.08000 + 
-0.05000 ^ 
0.01000 + 
0.23324 - 
0.40817-0 
0.34000 + 
0.22978 - 
0.40212 - 
0.33000 + 
-0.07937- 
-0.47624 - 
0.63000 + 
-0.16733 - 
-0.41833 - 
0.70000 + 
0.11314- 
-0.04243 ^ 
0.02000 + 



O.OlOOOi \ 

- 0.03000i 
O.OOOOOi / 

oooooA 

05831i 

00000 jy 

05745A 
05745i 
OOOOOi/ 

- o.oooooA 

- O.OOOOOi 
O.OOOOOi / 

- 0.08367i\ 

- O.OOOOOi 
O.OOOOOi J 
0.04243i \ 

- 0.056571 
O.OOOOOi J 



6.5 Discussion 



We obtained the similar result as that of Osawa and Nagaoka. Since we have not 
done any error estimation, we can not know how correct the result is. However, 
it is still meaningful because we achieved a good approximation for the result of 
Osawa and Nagaoka, and showed their result is likely to be the global optimum. 

In our experiment, the number of quantum states which appears on the surface 
of the optimum ball is 6, which is bigger than the number 3 shown in the paper by 
Osawa and Nagaoka. Although it is probably because of an error of the approxi- 
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mation, another possibility is that for the smaller number of points, the optimum 
is not correctly calculated. Actually we observed that smaller number of points 
as a constraint makes the optimization more unstable. It is because with fewer 
constraints, the dimension of the space where the variables can go becomes higher. 
Anyway, further observation is needed to know the reason for the difference. 

To make the approximation better, we have to make the step value A smaller. 
The problem is the computational complexity which rapidly increases as A gets 
smaller. For three level system, when A is 1/n, the complexity becomes O(n^). 

To avoid the computationally difficult case, our experiment was under the as- 
sumption that the optimum ball is determined by less than d? points. To remove 
the restriction is naturally considered as a next work. Although the direct compu- 
tation of the ball which pass through d? is to solve a d? — 1 dimensional non-linear 
equation, the dual problem by Legendre transform becomes a linear equation. Us- 
ing the dual problem is a probable option. 

Another problem is how much likely this method can be used for an arbitrary 
d-level system. For a d-level system, the complexity becomes 0{v?'^~'^), where 
A = 1/n. Because of lack of error estimation, we cannot say any concrete criterion 
for n, but it is likely that the computation for c? = 4 is far more difficult than the 
case d = 2>. 

As is mentioned by Welzl [121j . the possible way of improvement is to care for 
the order of the points. The original algorithm to compute the smallest enclosing 
ball is a Las Vegas algorithm and its computational complexity is based on the 
randomness of the order of points. Some heuristics might improve the actual per- 
formance of the computation, even if it does not improve the complexity expressed 
as an order. In particular, since the optimization process for fixed bounding points 
is very heavy in the proposed algorithm, the heuristics in the order of the points 
is more likely to contribute to the improvement than the usual Euclidean smallest 
enclosing ball problem. 
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6.6 Summary of this chapter 

We proposed a new algorithm to compute the Holevo capacity of a quantum channel 
which uses arranged version of Welzl's algorithm to compute the smallest enclosing 
ball problem. The proposed algorithm is for an arbitrary d-level system, and it is 
natural extension of the existing algorithm for one-qubit system. 

Although the proposed algorithm includes a non-linear optimization in its pro- 
cess which is likely to be very unstable, we have shown it is actually works for some 
examples of channels in 3-level system. We calculated the same examples as the 
ones by Osawa and Nagaoka [92], and certified the results by them are real optima 
though there had been a concern that those might be only local optima. The main 
merit of our algorithm is that it searches an optimum globally. 

Although proposed algorithm can practically approximate the result, analysis 
for its error is not done yet, and the error is not one-sided. Estimating the error 
and computing the lower-bound for the real value by the computed value are to be 
done. 
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Chapter 7 



Conclusion 

In this dissertation, we have given another geometric interpretation to a quantum 
state space. It indicates a connection among some quantum-related research fields. 

More precisely, what we have shown in this dissertation mainly consists of the 
following two parts: 

• some Voronoi diagrams in quantum state space coincides 

• Welzl's algorithm for the smallest enclosing ball problem is also applicable 
for a quantum state space. 

In this chapter, we summarize again the main results and discuss the meaning 
of them and future potential of development of related researches. 

7.1 Coincidences of Voronoi diagrams 

We proved that in one-qubit system, the Voronoi diagrams with respect to the 
divergence and Euclidean distance are the same. We also proved that in an n-level 
system for n > 3, that coincidence does not occur. 

In three or higher level system, restricting on a certain subspace, wc have in- 
vestigated some distances. We showed that with some unnatural parameterization, 
Voronoi diagrams with respect to the divergence and Euclidean distance coincide 
on the supspace. We also showed that Voronoi diagrams with respect to the diver- 
gence and the Fubini-Study distance coincide. 
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The existing only application of these fact is the algorithm by Hayashi et al. [U] 
to compute the Holevo capacity of a quantum channel. Nevertheless, we believe 
those result have clarified some aspect of the structure of a quantum state space. 

Our result shows thatthere is a connection among different distances which are 
used in different cotexts. Especially a significant thing is a connection between 
the distance used in quantum computation and the pseudo-distance which is used 
for quantum information theory. We regard we indicated a bridge between those 
research fields. 

From the viewpoint of computational geometry, what we have done can be a 
methodological hint. Generally in a situation that some measures are associated to 
a set, analyzing those Voronoi diagrams will some indication about the relation of 
those measures. Especially for a computational purpose, since a Voronoi diagram 
is a popular tool to approximate a continuous geometric object in a computer, 
comparing Voronoi diagrams is a reasonable start point for the discussion about 
how we can deal with a given continuous object with some associated metrics. 

7.2 Numerical computation of Holevo capacity 

We have shown that Welzl's algorithm for the smallest enclosing problem is also 
applicable to a quantum state space. Using Welzl's algorithm, we proposed a new 
algorithm to compute the Holevo capacity of a quantum channel. Although it leads 
to no conclusion about its complexity, we showed it is practical by an experiment. 

Although the experiment is for only few samples and under some restrictions, 
we showed the proposed algorithm really works for a real computation. We have 
not compared its practical performance to the existing algorithm by Osawa and 
Nagaoka [92], but it has at least one merit: there is no worry for local optima. 

Since our algorithm is an approximation by plotting points on a continuous 
geometric object, its performance is trade-off with its error. Consequently, for the 
future research, the improvement of "uniformness" of plotted points and the error 
estimation should be considered as a set. The proposed method to generated a 
"reasonably" uniform points as pure states is only intended to be implemented 
easily, and no mathematical analysis is given for its uniformness. To improve the 
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algorithm of that part, consideration about what error bound can be attained wih 
be necessary. 

Welzl's algorithm is classified as a Las Vegas algorithm, and its expected time 
of computation is linear. As is mention by Welzl himself |121j . the expected time is 
based on the randomness of the order of points, but some heuristics about ordering 
of points might improve the performance. To seek for heuristics specific for a 
quantum state space is one of the possible extensions of this research. 
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